AFRL-SN-RS-TR-2000-30  Vol  I  (of  two) 
Final  Technical  Report 
March  2001 


INCREASING  SECONDARY  DATA  OR  ENHANCED 
ADAPTIVE  DISPLACED  PHASE  CENTERED 
ANTENNA  APPROACHES 


Capraro  Technologies,  Inc. 
Rick  S.  Blum 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


20010507  067 


AIR  FORCE  RESEARCH  LABORATORY 
SENSORS  DIRECTORATE 
ROME  RESEARCH  SITE 
ROME,  NEW  YORK 


Although  this  report  references  (*)  limited  document (s)  on  pages  35  and  54 
no  limited  information  has  been  extracted* 

This  report  has  been  reviewed  by  the  Air  Force  Research  Laboratory,  Information 
Directorate,  Public  Affairs  Office  (IFOIPA)  and  is  releasable  to  the  National  Technical 
Information  Service  (NTIS).  At  NTIS  it  will  be  releasable  to  the  general  public, 
including  foreign  nations. 


AFRL-SN-RS-TR-2001-30  Vol  I  (of  two)  has  been  reviewed  and  is  approved  for 
publication. 


APPROVED:  ^ 

DAVID  B.  BUNKER 
Project  Engineer 


FOR  THE  DIRECTOR: 


ROBERT  G.  POLCE,  Chief 
Rome  Operations  Office 
Sensors  Directorate 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  Air  Force  Research 
Laboratory  Rome  Research  Site  mailing  list,  or  if  the  addressee  is  no  longer  employed  by 
your  organization,  please  notify  AFRL/SNRT,  26  Electronic  Pky,  Rome.  NY  13441- 
4514.  This  will  assist  us  in  maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a  specific 
document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing 
the  collection  of  mfonnation.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4301  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (OTOA-OISSh  Washington,  DC  20503. 


1.  (Leave  t/onk)  2.  REPORT  DATE 

MARCH  2001 


3.  REPORT  TYPE  AND  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

INCREASING  SECONDARY  DATA  OR  ENHANCED  ADAPTIVE  DISPLACED 
PHASE  CENTERED  ANTENNA  APPROACHES 


6.  AUTHOR(S) 
Rick  S.  Blum 


Final  Mar  97  -  Feb  98 


5.  FUNDING  NUMBERS 

C  -  F30602-97-C-0065 

PE  -  62702F 

PR  -  4506 

TA-  11 

WU-PJ 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Prime:  Capraro  Technologies,  Inc.  Sub:  Leigh  University 
311  Turner  Street  -  Suite  410  EECS  Department 
Utica  NY  13501  19  Memorial  Drive  West 

Bethlehem  PA  18015 


.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Air  Force  Research  Laboratory/SNRT 
26  Electronic  Pky 
Rome  NY  13441-4514 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

AFRL-SN-RS-TR-2001-30 
Vol  I  (of  two) 


11.  SUPPLEMENTARY  NOTES 

Air  Force  Research  Laboratory  Project  Engineer:  David  B.  Bunker/SNRT/(315)  330-2345 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


13.  ABSTRACT  (Maximum  200  words) 

A  Study  of  space-time  adaptive  processing  in  nonhomogeneous  environments  is  presented.  The  results  produced  are  logically 
organized  as  investigations  into  three  different  but  related  topics.  The  first  topic  can  be  summarized  as  an  investigation  of 
the  performance  of  partially  adaptive  space-time  adaptive  processing  schemes  using  measured  (MCARM)  data.  The  second 
topic  considers  the  development  of  a  new  set  of  space-time  processing  algorithms,  which  use  prior  knowledge  to  reduce  the 
number  of  parameters,  which  need  to  be  adapted.  The  new  schemes  were  tested  using  MCARM  data.  The  final  topic 
provides  analytical  realistic  cases  where  the  reference  or  training  data  may  have  different  statistics  from  the  data  in  the  ceil 
under  test.  Such  equations  have  not  yet  appeared  in  the  literature  and  they  are  useful  for  understanding  the  performance  of 
space-time  adaptive  processing  (STAP)  schemes  when  they  are  tested  with  MCARM  data. 


14.  SUBJECT  TERMS 


Space  Time  Adaptive  Processing  (STAP),  Adaptive  Array  Processing  (ASAP),  Estimation 
Using  Mismatched  Statistics,  Airborne  Radar  Signal  Processing 


15.  NUMBER  OF  PAGES 

120 

16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIHED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2-89  (EG) 

PrattfibfldbyANSI  Std.  239.18 

Dwignad  using  P^onn  Pro,  WHS/OIOR,  Oct  94 


Table  of  Contents 


Introduction 

Chapter  1  3 

1 .  Introduction  4 

2.  Some  Reduced  Complexity  STAP  Schemes  6 

2.1  General  STAP  Approach  7 

2.2  ADPCA  and  its  beamspace  version  8 

2.3  FTS  and  subarrayed  FTS  1 1 

2.4  EFA  and  subarrayed  EFA  11 

2.5  Joint-Domain  Localized  approach  (JDL)  12 

3.  Measured  Data  Results  13 

4.  Conclusions  13 

5.  References  35 

Chapter  2  38 

1 .  Introduction  39 

2.  General  approach  for  modified  scheme  40 

3.  Parameter  choice  45 

4.  Performance  comparison  48 

5.  Conclusion  33 

References  34 

Chapter  3  72 

Introduction  73 

Distribution  of  the  Test  Statistic  74 

Probability  of  detection  and  false  alarm  85 

Effects  of  Mismatch:  Numerical  results  87 

Extension  to  a  more  general  class  of  tests  89 

Steering  and  Covariance  Mismatch  90 

Steering  and  Covariance  Mismatch:  An  Explanation  for  the  good  performance  of 
ADPCA  97 

Conclusion  99 

References  99 

Appendix  191 


1 


Chapter  1 


List  of  Figures 


Figure  1)  ADPCA  sub-CPI  formation  9 

Figure  2)  Power  spectrum  plot  of  range  150  17 

Figure  3)  Performance  comparison  when  target  is  inserted  at  range  1 50, 

case  a.  18 

Figure  4)  Performance  comparison  when  target  is  inserted  at  range  150, 

caseb.  19 

Figure  5)  Performance  comparison  when  target  is  inserted  at  range  150, 

case  c.  20 

Figure  6)  Performance  comparison  when  target  is  inserted  at  range  150, 

case  d.  21 

Figure  7)  Performance  comparison  when  target  is  inserted  at  range  150, 

case  e.  22 

Figure  8)  Performance  comparison  when  target  is  inserted  at  range  150, 

case  f.  23 

Figure  9)  Performance  comparison  when  target  is  inserted  at  range  150, 

case  g.  24 

Figure  10)  Power  spectrum  plot  of  range  350  25 

Figure  11)  Performance  comparison  when  target  is  inserted  at  range  350, 

case  a.  26 

Figure  12)  Performance  comparison  when  target  is  inserted  at  range  350, 

case  b.  27 

Figure  13)  Performance  comparison  when  target  is  inserted  at  range  350, 

case  c.  28 

Figure  14)  Performance  comparison  when  target  is  inserted  at  range  350, 

case  d.  29 

Figure  15)  Performance  comparison  when  target  is  inserted  at  range  350, 

case  e.  30 

Figure  16)  Performance  comparison  when  target  is  inserted  at  range  350, 

case  f.  31 

Figure  17)  Performance  comparison  when  target  is  inserted  at  range  150, 

case  g.  32 

Figure  18)  Power  spectrum  plot  of  range  415  33 

Figure  19)  Power  spectrum  plot  of  range  500  34 


Chapter  2 

Figure  1)  Performance  comparison  for  target  inserted  in  different  range  bins 
and  with  several  different  normalized  Doppler  frequencies  for 
Modified  ADPCA  when  csfs^d  -  0.1,fi?=l,...,5.  55 

Figure  2)  Performance  comparison  for  target  inserted  in  different  range  bins 


11 


Figure  3) 

Figure  4) 

Figure  5) 
Figure  6) 

Figure  7) 

Figure  8) 

Figure  9) 

Figure  10) 

Figure  11) 


and  with  several  different  normalized  Doppler  frequencies ...  56 

Performance  for  modified  SMI  when  target  is  inserted  in  a 
particular  range  hin  with  a  particular  normalized  Doppler  frequency 5  7 
Performance  for  modified  SMI  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  frequency  5  8 
Average  energy  as  a  function  of  range  bin  59 

Performance  for  modified  ADPCA  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  ffequencydl 
Performance  for  modified  ADPCA  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  ffequency63 
Performance  for  modified  EFA  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  frequency65 
Performance  for  modified  EFA  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  frequency67 
Performance  for  modified  JDL  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  ffequency69 
Performance  for  modified  JDL  when  target  is  inserted  in  a 
particular  range  bin  with  a  particular  normalized  Doppler  frequency 71 


Chapter  3 


Figure  1 
Figure  2 
Figure  3 
Figure  4 
Figure  5 

Figure  6 
Figure  7 
Figure  8 


Probability  of  detection  variation  with  Var  {d} 

Probability  of  detection  variation  with  C-[(|)i,  ^2,  <1)3]- 
Probability  of  false  alarm  variation  with  Var  {d}. 

Probability  of  false  alarm  variation  with  C. 

Power  spectral  density  of  a  simple  model  [17]  for  ground  clutter 
and  mismatch  (extra  hump). 

Performance  variation  for  an  airborne  radar  example 
Probability  of  detection  variation  with  steering  mismatch  (offset). 
Probability  of  detection  variation  with  steering  mismatch  (not 
offset) 


103 

104 

105 

106 

107 

108 
109 

109 


List  of  Tables 


Chapter  1 

Table  1)  Test  cases  for  each  example 

Table  2)  The  three  best  schemes  for  all  cases  when  the  target  is  inserted  at 
range  150  (D  is  the  approximate  difference  between  the 
normalized  Test  statistic  at  some  other  range). 

Table  3)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted 
at  range  350. 

Table  4)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted 
at  range  415. 

Table  5)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted 
at  range  500. 

Table  6)  Parameters  for  comparison  test  in  this  study 

Chapter  3 

Table  1)  Parameters  of  the  power  spectral  density  model  in  Fig. 


INTRODUCTION 


This  report  documents  the  work  done  by  Lehigh  University  under  contract  F30602-97-C- 
0065-01.  This  work  involved  completing  tasks  4.1.1  and  4.1.2  of  the  original  contract. 
The  results  produced  are  logically  organized  as  investigations  into  three  different  but 
related  topics.  Each  of  these  topics  forms  one  chapter  in  this  report.  Each  of  these 
chapters  is  self-contained  and  includes  an  introduction  and  a  conclusion.  The  first 
chapter  includes  an  expanded  version  of  work  proposed  under  task  4.1.1.  This  work  can 
be  summarized  as  an  investigation  of  the  performance  of  partially  adaptive  space-time 
adaptive  processing  schemes  using  MCARM  data.  The  second  chapter  includes  an 
expanded  version  of  work  proposed  under  task  4.1.2.  This  work  can  be  summarized  as 
the  development  of  a  new  set  of  space-time  processing  algorithms  which  use  prior 
knowledge  to  reduce  the  number  of  parameters  which  need  to  be  adapted.  The  new 
schemes  were  tested  using  MCARM  data.  The  final  chapter  provides  analytical 
expressions  for  evaluating  the  performance  of  space-time  processing  schemes  in  realistic 
cases  where  the  reference  or  training  data  may  have  different  statistics  from  the  data  in 
the  cell-under-test.  Such  equations  have  not  yet  appeared  in  the  literature  and  they  are 
useful  for  understanding  the  performance  of  space-time  adaptive  processing  schemes 
when  they  are  tested  with  MCARM  data.  The  final  chapter  has  results  which  are 
applicable  to  both  task  4.1.1  and  task  4. 1 .2. 
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Chapter  1  :  Task  4.1.1 


1.  INTRODUCTION 

In  airborne  radar,  target  detection  is  often  limited  by  ground  clutter  and  other 
forms  of  interference.  Platform  motion  causes  Doppler  shifts  in  interference  returns, 
making  Doppler  filtering  alone  ineffective.  In  such  cases  Space-Time  Adaptive 
Processing  (STAP)  offers  a  potential  solution. 

STAP  has  been  an  active  research  topic  for  at  least  the  last  two  decades.  Much  of 
the  interest  was  generated  by  the  research  in  [1]  and  [2].  Since  then  several  algorithms 
have  been  proposed  and  evaluated  using  simulated  radar  data.  With  the  recent 
improvements  in  phased  array  antenna  and  digital  signal  processing  technology,  a  STAP- 
based  airborne  radar  system  is  becoming  an  attractive  alternative  for  detection  of  airborne 
targets  competing  with  strong  clutter  returns  and  jamming,  as  compared  to  classical  low- 
sidelobe  beamforming  [3]. 

Current  STAP  research  efforts  [4]  are  focused  on  a  number  of  interesting  issues. 
Performance  evaluation  of  several  competing  STAP  approaches  is  the  topic  of  this  report. 
In  most  previous  research,  STAP  algorithms  were  evaluated  using  optimistic  simulated 
data  or  by  manipulating  stationary  platform  measurements  to  emulate  motion.  While 
simulated  data  is  very  useful  in  the  development  and  analysis  of  algorithms,  a  more 
complete  evaluation  includes  using  actual  recorded  radar  data.  Thus,  in  this  report,  we 
compare  various  STAP  schemes  using  measured  airborne  data.  A  general  formulation  of 
a  useful  class  of  STAP  processing  approaches  is  developed.  This  formulation,  which 
includes  most  linear  processing  schemes,  unifies  the  techniques  considered. 

Assume  the  radar  transmits  a  coherent  burst  of  M  pulses  at  a  constant  pulse 
repetition  frequency  fr=  1/Tr,  where  Tr  is  the  pulse-repetition-interval  (PRI).  The  time 
interval  over  which  the  waveform  returns  are  collected  is  commonly  referred  to  as  the 
coherent-processing-interval  (CPI).  The  radar  antenna  used  is  an  array  antenna  with  N 
identical  elements.  Denote  the  observation  corresponding  to  the  i***  antenna  element,  the 
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pulse  and  the  range  cell  as  Xy,k.  It  is  convenient  to  denote  the  data  associated  with 
the  k***  range  cell  as 


(1) 


where  denotes  the  transpose  of  the  vector  a.  We  will  refer  to  Xk  as  a  space-time 
snapshot. 

A  simple  set  of  STAP  schemes  that  have  been  suggested  can  be  represented  as 
an  inner  product  of  the  conjugate  of  a  weight  vector  w  and  the  vector  Xk-  This  inner 
product 

Z  =  W%  (2) 

produces  the  complex  quantity  z  whose  magnitude  is  often  compared  to  a  threshold  to 
make  a  decision  regarding  target  presence  or  absence.  The  weight  vector  w  may  depend 
on  the  estimated  interference-plus-noise  environment  and  on  the  target  of  interest.  In  the 
well-know  sample  matrix  inversion  (SMI)  algorithm  [2]  for  example,  a  popular  fully 
adaptive  algorithm,  the  weight  vector  is  given,  to  within  a  scale  factor,  as 

w  =  R'V  ® 

where  R  is  the  estimated  interference-plus-noise  covariance  matrix.  The  estimate  is 
based  on  a  set  of  reference  data,  typically  chosen  from  the  surrounding  range  cells.  V  is 
the  normalized  target  response  (Xk  for  a  target  observed  without  clutter  or  noise)  [5].  In 
the  case  where  the  interference  statistics  are  known  or  the  estimated  covariance  matrix  is 
exactly  equal  to  the  true  covariance  matrix,  SMI  can  achieve  optimal  performance.  A 
fully  adaptive  STAP  scheme,  like  SMI,  requires  the  formation  of  an  NM  by  NM 
covariance  matrix.  For  moderate  M  and  N,  the  computational  cost  of  the  computation  of 
R*^  becomes  excessive  for  real-time  implementation.  Further,  schemes  which  estimate  R 
typically  require  a  large  set  of  independent  and  identically  distributed  (iid)  reference  data 
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vectors  to  achieve  an  accurate  estimate  for  moderate  M  and  N.  This  requirement  may  be 
unrealistic,  since  measurements  [6]  indicate  that  multi-channel  airborne  radar  clutter  data 
is  often  severely  non-homogeneous.  For  this  reason  the  reference  data  set  available  for 
estimation  of  clutter  statistics  is  usually  quite  small.  As  a  result,  reduced  complexity 
approaches  called  partially  adaptive  STAP  (PASTAP)  have  been  developed  whose 
computational  cost  and  reference  data  requirements  are  more  realistic.  Some  examples  of 
reduced  complexity  approaches  are  given  in  the  next  section. 

It  is  important  to  know  how  different  PASTAP  algorithms  perform  for  realistic  cases. 
This  issue  is  studied  in  this  report.  In  Section  2  we  define  a  general  PASTAP  scheme  and 
give  a  detailed  description  of  several  specific  approaches.  Performance  comparisons 
based  on  measured  airborne  radar  data  are  presented  in  Section  3.  Conclusions  are  given 
in  Section  4. 

2.  Some  Reduced  Complexity  STAP  Schemes 


STAP  is  an  active  research  area  and  new  schemes  are  continually  being  developed. 
In  order  to  compare  schemes,  a  standard  terminology  is  useful.  Here,  we  will  mainly 
follow  the  terminology  used  in  [5].  We  caution  the  reader  that  other  terminology  also 
appears  in  the  literature.  We  first  define  a  general  formulation  of  a  PASTAP  approach 
which  encompasses  many  existing  algorithms.  Next  we  describe  eight  specific 
approaches  which  are  included  in  the  general  formulation.  They  are 

1.  Adaptive  Displaced  Phase-Centered  Antenna  (ADPCA) 

2.  subarrayed  ADPCA  (BDPCA) 

3.  Beamspace  ADPCA  (BeamAD) 

4.  Factored  Post  Doppler  (FTS) 

5.  subarrayed  FTS  (BFTS) 

6.  Extended  Factored  Approached  (EFA) 

7.  subarrayed  EFA  (BEFA) 
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8.  Joint-Domain  Localized  Approach  (JDL) 


2.1  General  STAP  Approach 


Consider  the  transformations 


Xk(p)=  (Ap  0  Bp]^X|j ;  ps:0,l,2,...,P-l  (4) 

where  Xr  is  the  space-time  snapshot  from  the  k”*  range  cell  and  Ap  and  Bp  are  scheme- 
dependent  matrices.  The  operations  in  (4)  can  be  interpreted  as  a  pre-processor  applied 
to  the  received  signals.  This  pre-processing  generates  data  for  the  adaptive  processing  to 
follow.  Note  that  P  vectors  are  produced  by  the  operations  in  (4).  Typically,  the  pre¬ 
processing  in  (4)  performs  a  coordinate  transformation  and  a  selection  operation. 

We  describe  the  adaptive  processing  on  the  p***  vector  as 

y,(p)=S"R;’(p)X,(p)A>  (5) 

where 


1  k+Q/2+1  -  -  ,  vu 

R|.(p)  =  K  I  X,(p)x,(p)" 

Q  i=k-Q/2-l,i;‘k-IJe,k+l 


(6) 


and  S  is  a  scheme-dependent  steering  vector.  Rk(p)  is  the  interference-plus-noise 
covariance  matrix  estimated  from  Q  adjacent  range  cells,  excluding  the  cell-under-test 
and  the  two  closest  range  cells.  The  estimate  in  (6)  is  the  maximum  likelihood  estimate  if 
the  data  samples  from  the  Q  adjacent  range  cells  are  independent  and  identically 
distributed  with  the  same  complex  Gaussian  distribution  as  the  interference-plus-noise  in 
the  cell-under-test.  It  is  straight-forward  to  generalize  (6)  to  allow  the  range  cells  to  be 
selected  differently,  including  selections  which  are  non-symmetric  about  the  cell-under¬ 
test.  We  do  not  discuss  these  generalizations  here  since  we  use  only  the  approach 
suggested  in  (6)  in  our  tests.  The  term  O  in  the  denominator  of  (5)  is  the  normalization  to 


provide  CFAR  in  homogeneous  clutter  and  is  given  by 


4>  =  Vs”r;;'s 


(7) 
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In  different  schemes,  yk(p)  may  or  may  not  be  the  final  output  of  interest.  If  yk(p) 
is  the  final  output  of  interest,  its  magnitude  will  be  compared  to  a  threshold  to  decide  if 
target  is  present.  For  cases  where  y^Cp)  will  be  processed  further,  we  assemble  the 
complex  outputs  from  each  adaptive  processor  as 

Yk  =  [y  K  (0),  y ,  (1), . . . ,  y,  (P  -  l)f  (8) 

and  compute 

Zk.m=f”Y,  (9) 

which  we  call  post-processing  (after  adaptive  processing).  Typically,  fn,  is  the  m*'* 
column  of  a  PxP  filter  matrix  F,  and  Zk,™  is  the  final  output  whose  magnitude  will  be 
compared  to  a  threshold  to  produce  a  decision. 


2.2  ADPCA  and  its  beamspace  version 

Define  a  set  of  P  sub-CPIs  Xk(p),  p=0,...,P-l  in  the  k***  snapshot.  Each  sub-CPI 
contains  possible  target  returns  from  Kj  pulses  and  all  N  elements.  Fig.  1  shows  two 
different  ways  to  form  the  sub-CPIs.  As  indicated  in  Fig.  1,  implementation  (a)  does  not 
overlap  pulses.  Given  M  pulses  in  a  CPI  where  M  can  be  divided  by  Kt,  implementation 
(a)  generates  P=M/Kt  sub-CPIs.  The  0“*  sub-CPI  consists  of  pulses  0,  ...,  Kfl  and  the 
p"*  sub-CPI  consists  of  pulses  pKt, pKt+Kfl.  Implementation  (b)  forms  the  sub-CPIs 
by  using  the  same  pulse  returns  in  several  sub-CPIs.  Given  M  pulses  in  a  CPI, 
implementation  (b)  generates  P=:M-K,+1  sub-CPIs.  The  0“'  sub-CPI  consists  of  pulses  0, 
Kfl  and  the  p“*  sub-CPI  consists  of  pulses  p,  p+K,-!.  In  Fig.l,  Kt  is  set  to  3 
and  in  implementation  (b),  neighboring  sub-CPIs  overlap  2  pulses.  Of  course,  other 
overlaping  pulse  configurations  are  possible. 

The  pre-processing  we  have  just  described  can  be  put  into  the  framework  of  (4). 
Bp  is  set  to  In  which  is  an  NxN  identity  matrix  and 
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(10) 


”p(K,-h><K, 

'k. 

^(M-K,-pK,+ph><K, 


where  the  notation  Oqxm  refers  to  an  qxm  matrix  of  zeros,  h  indicates  the  number  of 
pulses  which  are  overlapped.  In  implementation  (a)  h  is  set  to  be  zero  and  in 
implementation  (b)  h  is  set  to  be  Kfl.  Ap  is  an  MxKt  selection  matrix. 
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Figure  1)  ADPCA  sub-CPI  formation 


The  adaptive  processing  in  ADPCA  is  described  by  (5)  with  the  steering  vector 


S  =  S,  ®  Sg 


(11) 


where  S*  is  the  Nxl  spatial  steering  vector  [5]  and  St  is  a  KtXl  vector,  which  is  composed 
of  the  binomial  coefficients,  with  each  coefficient  altered  in  sign  (start  with  positive).  As 
a  particular  example,  we  have 


S,  =  (l,-2,iy 


(12) 
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for  a  three  pulse  case. 

Typically,  post-processing  as  described  in  (8)  and  (9)  is  employed  in  ADPCA.  In 
ADPCA  F  is  a  matrix  corresponding  to  a  Doppler  filter  bank,  and  is  the  filter 
corresponding  to  the  m“’  Doppler  frequency.  A  DFT  matrix  implements  the  Doppler  filter 
bank.  Thus,  a  fast  fourier  transform  efficiently  computes  (9)  for  all  m.  IZk^l,  the 
magnitude  of  the  final  output  for  Doppler  bin  m,  is  compared  to  a  threshold  to  test  for  a 
target  in  Doppler  bin  m. 

Subarrayed  ADPCA  and  beamspace  ADPCA  are  both  beamspace  versions  of 
ADPCA.  Bp  is  used  as  a  beamformer  matrix  to  produce  K*  <  N  beams.  We  focus  on 
Ks=3  in  this  report. 

In  subarrayed  ADPCA  the  beamformer  matrix 


is  employed  in  (4) ,  where  g  =  (go,  gi,  g2, . . gN-.i)'^  is  an  N'xl  vector  and  N'=  N-Ks+1. 
The  vector  g  can  be  any  of  the  popular  window  functions  from  the  DSP  or  radar  literature 
[7,8].  Only  results  which  use  the  uniform  window  are  provided  here.  In  our  preliminary 
tests,  the  uniform  window  appeared  to  outperform  the  other  windows  we  tried  so  we 
focused  on  this  choice. 

Each  vector  produced  by  (4)  will  be  adaptively  processed  as  in  (5)  where 
S  =  (lKt®G)”S,  (14) 

with  Se  chosen  as  the  steering  vector  defined  by  (1 1).  The  post  processing  is  the  same  as 
in  ADPCA. 
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Beamspace  ADPCA  is  identical  to  subarrayed  ADPCA,  but  with  Bp=  G  =[  fN,i, 

•  •  •»  fN,Ks  ].  where  Tnj  j=l>  •  •  •,  Ks  are  those  columns  of  NxN  DPT  matrix  corresponding 
to  Ks  particular  angle  bins. 


2.3  FTS  and  subarrayed  FTS 

In  factored  post-Doppler  STAP  [5],  Doppler  processing  is  first  performed  on  each 
spatial  channel.  Let  the  Doppler  filter  applied  to  the  p“*  spatial  channel  be  fp  and  assume 
fp  is  chosen  from  the  set  of  possible  Doppler  filters  in  the  MxM  matrix  Fm=[  fo>  fij  -..jfivi- 
i].  Then  the  pre-processing  is  described  by  (4)  with  Ap=fp,  P=M,  and  Bp=lN.  This  pre¬ 
processing  transforms  the  signal  into  Doppler  space.  In  this  case  p  indicates  the  index  of 
Doppler  bin  in  question.  Next,  the  adaptive  processing  in  (5)  is  employed  with  the 
steering  vector  defined  as  in  (1 1)  with  Kt  =1.  Here  S=Ss,  since  Kt  is  1.  As  for  most  of  the 
STAP  schemes  we  discuss,  tapering  could  be  applied  to  the  steering  vector  [5].  Post¬ 
processing  is  not  usually  employed  and  lyk(p)l  is  compared  to  a  threshold  to  test  for  a 
target  in  the  p***  Doppler  bin. 

Subarrayed  FTS  is  identical  to  FTS  but  with  Bp=  G  as  described  in  (13).  The 
adaptive  processing  is  as  in  FTS  with  the  steering  vector  defined  in  (14).  Here  Se  is  the 
steering  vector  that  is  used  in  the  adaptive  processing  of  FTS.  Post  processing  is  not 
employed. 


2.4  EFA  andsubarrayed  EFA 

The  extended  factored  approach  (EFA)  [9]  is  a  slight  extension  of  the  factored 
post-Doppler  approach.  In  EFA,  adaptive  processing  is  applied  to  several  adjacent 
Doppler  bins  instead  of  just  one.  Thus,  the  pre-processing  performs  both  transformation 
and  selection.  In  the  case  considered  here,  where  the  scheme  adapts  over  3  adjacent  bins, 
the  pre-processing  can  be  described  as  in  (4)  with  Ap  =  Jp  =  [  fp.i,  fp,  fp+i]  and  Bp=lN. 


The  other  quantities  are  set  in  a  similar  manner  as  for  factored  post-Doppler  STAP 
but  with  the  new  Ap.  For  example,  the  steering  vector  is  obtained  from 

S  =  (Ap®Bp)fV  (15) 

where  V  is  the  normalized  target  response  as  used  in  (3).  For  example,  if  fp  correspond  to 
the  target  Doppler  frequency  under  consideration,  the  steering  vector  can  be  defined  as  in 
(11)  with  St  =  [0,1,  Of  . 

In  subarrayed  EFA,  the  pre-processing  is  the  same  as  with  EFA  but  with  Bp  set  as 
in  (13).  The  adaptive  processing  is  defined  as  in  (5)  with  S  as  in  (15). 

2.5  Joint-Domaint  Localized  approach  (JDL) 


In  JDL  [10],  the  pre-processor  performs  two  dimensional  transformation  and 
selection.  The  data  is  transformed  from  the  space-time  domain  into  the  angle-Doppler 
domain.  This  pre-processing  can  be  described  as  in  (4),  with  Ap=[fm,i,  ...,  fni,Kt], 
where  f„j  j=l,  ...,  Kt  are  Kt  columns  of  an  MxM  DFT  matrix  and  with  Bp=  G=[  f„4, 
fn^,  ...,  f„,Ks  ],  where  fnj  j=l,  ...,  Ks  are  Kg  columns  of  NxN  DFT  matrix.  As  for  FTS 
and  EFA,  only  the  post  processing  corresponding  to  a  single  p  must  be  calculated  to  test 
for  a  target  at  a  particular  normalized  Doppler  and  spatial  frequency. 

If  we  focus  on  Kt  =  Kg  =3  and  consider  the  case  where  the  target  to  be  detected 
has  the  Doppler  frequency  corresponding  to  fm,2  and  the  spatial  frequency  corresponding 
to  fn,2.  then  the  adaptive  processing  in  (5)  is  used  with 

S  =  [0  1  0]®[0  1  O]  (16) 

More  precisely,  the  steering  vector  has  all  its  entries  equal  to  zero  except  for  the 
one  corresponding  to  the  spatial  and  Doppler  frequency  of  the  target.  No  post-processing 
is  employed  for  JDL. 
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3*  MEASURED  DATA  RESULTS 


To  test  the  STAP  algorithms  described  in  section  2,  we  use  data  from  the  Multi- 
Channel  Airborne  Radar  Measurements  (MCARM)  database,  flight  5,  acquisition  575. 
See  [11,12]  for  detailed  information  about  the  MCARM  program  and  the  data.  For  each 
experiment,  a  single  target  signal  with  amplitude  0.05,  a  particular  normalized  Doppler 
frequency  and  a  particular  spatial  frequency  was  inserted  in  a  particular  range  bin. 
Reference  data  are  selected  from  consecutive  range  cells  on  each  side  of  the  cell-under¬ 
test,  excluding  the  cell-under-test  and  the  two  closest  cells.  We  employ  normalized  test 
statistics  as  in  (5),  which  provide  a  constant  false  alarm  rate  (CFAR)  characteristic  for 
homogenous  clutter  [13].  The  parameters  used  in  each  of  the  STAP  algorithms  studied 
are  given  in  the  Appendix. 

For  each  example,  we  provide  plots  of  the  magnitude  of  the  normalized  test 
statistics  for  a  set  of  range  bins  including  the  target  range.  We  judge  a  scheme  by  how 
large  the  test  statistic  is  at  the  target  range  in  comparison  to  other  ranges. 

In  the  first  example,  we  inserted  a  target  at  range  bin  150  for  the  cases  shown  in 
Table  1,  The  location  of  the  targets  and  an  estimate  of  the  clutter  (plus  noise)  power 
spectral  density  (psd)  are  illustrated  in  Fig.  2.  As  visible  from  Fig.  2,  the  psd  estimate 
used  is  rather  crude  and  is  provided  to  give  a  rough  description  of  the  clutter 
environment.  In  the  estimate,  the  periodogram  method  [14]  is  used  with  blackman 
windowing.  The  operations  are  applied  just  on  the  cell-under- test  data  (no  neighboring 
range  cells  are  averaged).  We  found  that  some  of  the  artifacts  can  be  removed  by 
averaging,  but  this  was  not  considered  necessary  in  this  case. 
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Case 

Normalized  Doppler 

Frequency 

Spatial 

Frequency 

a 

-0.2656 

-0.2656 

b 

-0.0312 

-0.0312 

c 

0.125 

0.203 

d 

-0.0312 

-0.1875 

e 

-0.0312 

0.203 

f 

0.3593 

-0.1875 

g 

-0.1875 

0.3593 

Table  1)  Test  cases  for  each  example. 

A  summary  of  the  results  is  given  in  Table  2.  Fig.  3  through  Fig.  9  present  the 
results  for  most  of  the  schemes  tested.  FTS  and  subarrayed  FTS  generally  perform 
poorly,  so  their  results  are  not  shown.  In  this  example,  generally  JDL  provides  best 
results.  JDL  is  best  in  every  case  except  cases  d  and  e.  In  case  d,  BeamAD  slightly 
outperforms  JDL  and  ADPCA,  but  the  difference  is  quite  small.  In  case  e,  EFA 
outperforms  the  others,  however  JDL  also  performs  well.  The  apparently  large 
interference  indicated  by  Fig.  2  near  the  normalized  Doppler  frequency  -0.0312  could  be 
the  reason  why  JDL  is  not  best  in  these  two  cases.  These  are  apparently  difficult  cases 
where  no  scheme  can  really  excel.  Table  2  shows  that  case  b  is  also  quite  difficult,  but 
here  JDL  performs  much  better  than  the  other  schemes.  The  extra  clutter  ridges  in  Fig.  2 
are  discussed  in  [15]  and  [16]. 

In  the  second  example,  we  inserted  a  target  at  range  bin  350.  We  present  results 
for  the  same  cases  in  Table  1.  The  best  three  schemes  for  all  the  cases  are  given  in  Table 
3.  The  location  of  the  targets  and  an  estimate  of  the  clutter  psd  is  given  in  Fig.  10.  The 
normalized  test  statistics  for  the  six  best  schemes  tested  for  each  case  in  Table  1  are  given 
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in  Figures  1 1  through  17.  The  results  indie  site  thut  none  of  the  schemes  3lways 
outperforms  all  the  others.  However,  post-Doppler  algorithms  are  generally  better  than 
pre-Doppler  algorithms.  Either  JDL  or  EFA  were  best  in  all  but  case  b.  In  case  b,  where 
the  target  is  inserted  in  the  largest  clutter  of  all  cases,  BEFA  is  only  slightly  better  than 
EFA  and  JDL. 

Next,  we  inserted  a  target  at  range  bin  415.  We  present  results  for  the  same  cases 
as  in  previous  examples.  The  location  of  the  targets  and  an  estimate  of  the  clutter  psd  is 
given  in  Fig.  18.  A  summary  of  the  results  is  provided  in  Table  4.  Here,  for  cases  b,  c,  d, 
f,  and  g  either  JDL,  or  BEFA  provide  best  performance.  In  the  other  cases,  EFA  is  best 
and  JDL  also  performs  well.  JDL  performs  well  in  every  case  except  case  b.  Even  in  this 
case,  its  performance  is  adequate.  EFA  and  its  beamspace  version  are  best  in  some  cases 
and  near-best  in  others.  ADPCA  and  its  beamspace  version  give  good  performance  in 
number  of  cases,  but  these  schemes  were  never  best  in  this  example. 

Finally,  we  test  the  same  cases  as  before  when  the  target  is  inserted  in  range  bin 
500.  The  locations  of  the  targets  and  an  estimate  of  the  clutter  psd  is  given  in  Fig.  19.  A 
summary  of  the  results  is  given  in  Table  5.  Again  JDL,  EFA  and  BEFA  outperform  the 
other  algorithms  except  for  case  e  where  ADPCA  is  best. 

4.  CONCLUSIONS 

In  our  tests,  JDL  and  EFA  generally  perform  very  well.  Subarrayed  EFA  also 
shows  good  performance  in  many  situations.  The  common  element  these  schemes  share 
is  post  Doppler  processing.  This  type  of  processing,  used  in  the  correct  way,  appears  to 
be  superior  when  using  measured  data.  This  appears  to  be  the  major  result  of  this  study. 
When  there  is  a  strong  interference  near  the  target,  case  d  and  e  in  Table  5  for  example, 
ADPCA  performs  well  and  sometimes  outperforms  all  the  other  schemes.  The  reason 
appears  to  be  related  to  the  extra  whitening  provided  by  its  steering  vector. 

This  is  one  of  the  few  STAP  studies  we  have  seen  which  uses  real  airborne  radar 
measurements  and  thus  we  feel  these  results  are  interesting.  Judging  performance  using 


measured  data  is  difficult  since  one  can’t  directly  extract  probability  of  detection  and 
probability  of  false  alarm,  the  accepted  measure  of  performance  for  radar  signal  detection 
problems.  However,  measured  data  studies  are  still  important  to  obtain  a  more  complete 
assessment  of  performance.  In  the  current  study  we  present  comparisons  only  for  a 
particular  configuration  of  each  algorithm.  We  have  tried  to  pick  the  most  popular  or  at 
least  a  reasonable  configuration  for  each  algorithm.  The  results  may  be  different  for 
different  configurations.  Further  we  have  tested  the  algorithms  for  many  different  range 
bins  for  two  different  data  sets  (flights)  and  based  on  these  cases  the  results  given  here 
appear  to  be  representative.  However,  the  results  could  be  much  different  for  some  data 
sets  we  have  not  tested. 

We  have  obtained  results  for  a  fairly  large  number  of  different  cases.  Due  to  space 
limitations  we  have  provided  only  a  limited  set  of  these  in  this  report.  We  believe  that 
further  study  using  measured  data  is  needed  and  we  hope  to  see  more  papers  by  other 
research  groups  on  this  topic.  Hopefully  other  data  sources  will  become  available  and 
these  can  be  compared  to  results  obtained  using  the  MCARM  database.  Most 
importantly,  we  believe  that  techniques  for  assessing  performance  with  measured  data  is 
itself  a  topic  which  deserves  attention. 
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Normalized  doppler  frequency 


Normalized  spatial  frequency 


Figure  2)  Power  spectrum  plot  of  range  150 


Case 

Normalized  Doppler 

The  3  best 

H 

Figure 

Frequency 

■ESSSH 

schemes 

a 

JDL 

36 

-0.2656 

-0.2656 

BEFA 

25 

3 

BeamAD 

16 

-0.0312 

-0.0312 

JDL 

10 

4 

c 

JDL 

30 

0.125 

0.203 

EFA 

12 

5 

BeamAD 

8 

-0.1875 

JDL 

7 

6 

ADPCA 

6 

e 

EFA 

20 

-0.0312 

JDL 

12 

7 

BeamAD 

5 

f 

JDL 

43 

0.3593 

-0.1875 

EFA 

38 

8 

BeamAD 

30 

g 

JDL. 

38 

-0.1875 

0.3593 

BEFA 

18 

9 

BeamAD 

7 

Table  2)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted  at  range  150 
(D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target  and  the 
largest  peak  in  the  normalized  test  statistic  at  some  other  range). 
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AOPCA 


BOPCA 


a)  ADPCA 


b)  BDPCA 


c)  JDL 


d)  BeamAD 


e)  EFA 


f)  BEFA 


Figure  3)  Performance  comparison  when  target  is  inserted  at  range150,  case  a. 
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e)  EFA  f)  BEFA 


Figure  4)  Performance  comparison  when  target  is  inserted  at  range  150,  case  b. 
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AOPCA 


80PCA 


e)  EFA  f)  bEFA 


Figure  5)  Performance  comparison  when  target  is  inserted  at  range  150,  case  c. 
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e)  EFA  0  BEFA 

Figure  6)  Performance  comparison  when  target  is  inserted  at  range  150,  case  d. 
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e)  EFA  f)  befa 


Figure  7)  Performance  comparison  when  target  is  inserted  at  range  150,  case  e. 
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e)  EFA  f)  BEFA 


Figure  8)  PerfotTnance  comparison  when  target  is  inserted  at  range  150,  case  f. 
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e)  EFA 


Figure  9)  Performance  comparison  when  taro* 


2‘ 


Normalized  Doppler  frequency 


EFA 

JDL 

ADPCA 


EFA 

JDL 

ADPCA 


EFA 

JDL 

BeamAD 


JDL 

BEFA 

BeamAD 


Table  3)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted  at  range  350. 
(D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target  and  the 
largest  peak  in  the  normalized  test  statistic  at  some  other  range.). 
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a)  ADPCA 


b)  BDPCA 


c)  JDL 


d)  BeamAD 
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e)  EFA 


f)  BEFA 


Figure  11)  Performance  comparison  when  target  is  inserted  at  range  350,  case  a. 
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e)  EFA  f)  BEFA 

Figure12)  Performance  comparison  when  target  is  inserted  at  range  350,  case  b. 
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AOPCA 


BOPCA 


e)  EFA  f)  beFA 

Figure  13)  Performance  comparison  when  target  is  inserted  at  range  350,  case  c. 
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a)  ADPCA 


b)  BDPCA 


e)EFA 


f)BEFA 


Figure  15)  Performance  comparison  when  target  is  inserted  at  range  350,  case  e, 


a)  ADPCA 


b)  BDPCA 


NormaKzed  doppler  frequency 


•  T  J  \ 

-  0.5  0  0.5 

Normalized  spatial  frequency 


Figure  18)  Power  spectrum  plot  of  range  415 


Table  4)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted  at  range 
415.  (D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target 
and  the  largest  peak  in  the  normalized  test  statistic  at  some  other  range) 
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NormaGzed  spatial  frequency 


Figure  19)  Power  spectrum  plot  of  range  500 


Case 

Normalized  Doppler 
Frequency 

Spatial  Frequency 

The  3  best 
schemes 

D 

a 

JDL 

32 

-  0.2656 

-  0.2656 

EFA 

13 

ADPCA 

7 

JDL 

22 

-0.0312 

-0.0312 

EFA 

20 

BEFA 

15 

JDL 

36 

0.125 

BEFA 

12 

EFA,  BeamAD 

7 

d 

JDL 

26 

EFA 

20 

ADPCA 

10 

e 

ADPCA 

22 
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500.  (D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target 
peak  and  the  largest  peak  in  the  normalized  test  statistic  at  some  other  range) 
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APPENDIX 


Scheme 

M 

N 

Kt  1 

Ks  1 

Ap 

G 

DOF 

ADPCA 

128 

22 

■ 

[  0px3, 13,  0(i25.p)x3 

I22 

66 

EFA 

128 

22 

3 

■ 

[f  128, Dp-1  if  128,Dpif  128.Dp+l] 

I22 

66 

FTS 

128 

22 

1 

- 

tl28,Dp 

I22 

22 

Beamspace 

ADPCA 

128 

22 

3 

3 

[  0px3. 13,  0(125-p)x3 

See 

Note 

9 

Subarraying 

ADPCA 

128 

22 

3 

3 

[  0px3, 13,  0(125-p)x3  r 

*Q 

9 

Subarraying 

EFA 

128 

22 

3 

3 

[fl28,Dp-1ifl28,Dpifl28,Dp+l] 

*G 

9 

Subarraying 

FTS 

128 

22 

3 

3 

f  128, Dp 

*G 

3 

JDL 

128 

22 

3 

3 

[fl28,Dp-1ifl28,Dpifl28,Dp+l] 

See 

Note 

9 

Table  6)  Parameters  for  comparison  tests  in  this  study. 

DOF:  number  of  degrees  of  freedom.  It  is  equal  to  the  Q  in  (6)  in  this  paper. 

*  G  is  the  beamforming  matrix  which  is  formed  as  described  in  (13). 

Note:  MCARM  data  is  not  collected  by  a  uniformly  spaced  linear  antenna  array, 
so  its  beamforming  matrix  is  relatively  complicated.  This  G  matrix  used  in 
beamspace  ADPCA  and  JDL  can  be  obtained  by  using  measured  steering 
vector  information  provided  with  the  MCARM  database.  See  the  discussion  on 
modSAat[11]. 
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Chapter  2:  Task  4.1.2 
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1.  Introduction 


Adaptive  space-time  processing  schemes  for  airborne  radar  have  received  significant 
attention  recently  [1].  Typically  these  schemes  involve  using  observations  from 
neighboring  range  cells,  called  reference  data,  to  estimate  the  covariance  matrix  of  the 
clutter-plus-interference  in  the  range  cell-under-test.  The  estimated  covariance  matrix  is 
then  used  to  define  a  weight  vector  in  a  linear  processing  scheme.  The  weight  vector  is 
chosen  with  the  hope  that  it  will  suppress  the  clutter-plus-interference  in  the  range  cell- 
under-test.  If  the  estimated  covariance  matrix  accurately  represents  the  true  covariance 
matrix  for  the  cell-under-test  then  most  STAP  schemes  perform  well.  There  are  many 
reasons  why,  in  practice,  this  may  not  occur.  The  data  taken  from  surrounding  range  cells 
often  have  different  statistics  from  data  taken  from  the  cell-under-test.  Further,  the  amount 
of  reference  data  available  which  can  be  used  in  the  estimation  is  usually  quite  limited. 
Often  this  is  a  consequence  of  nonstationarity. 

To  overcome  these  difficulties,  reduced-dimension  STAP  schemes  [1,  2,  3,  4,  5]  have 
been  devised  which  have  shown  promise.  In  our  own  research  [3,  4,  5],  we  have  shown 
that  a  particular  approach  suggested  by  Rome  Laboratory  engineers,  the  Adaptive 
Displaced  Phase  Centered  Antenna  (ADPCA)  approach,  provides  good  performance  in 
some  cases  with  real  data.  Other  researchers  [1,4]  have  shown  that  the  Extended  Factored 
Approach  (EFA)  and  the  Joint-Domain  Localized  (JDL)  approach  can  also  provide  good 
performance  in  some  cases.  In  this  research,  a  new  methodology  for  constructing  robust 
STAP  schemes  is  proposed.  In  order  to  demonstrate  the  idea  we  use  a  simple  model  for  a 
case  with  ground  clutter  only  and  we  show  that  the  resulting  scheme  can  provide  good 
performance  in  tests  with  measured  radar  data. 
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In  section  2,  we  introduce  the  new  approach  and  describe  its  application  to  a  set  of  STAP 
algorithms.  In  section  3,  we  describe  how  to  choose  free  parameters  in  the  approach.  In 
section  4,  we  discuss  the  application  of  the  approach  to  four  specific  STAP  algorithms. 
Numerical  performance  results  are  given  for  a  measured  radar  data  set.  Conclusions  are 
given  in  section  5. 


2.  General  approach  for  modified  scheme 


Here  we  study  a  new  set  of  space-time  processing  algorithms  for  airborne  radar.  In  the 
spirit  of  ADPCA  and  other  reduced-dimension  algorithms,  our  new  algorithms  attempt  to 
estimate  fewer  parameters  than  fully-adaptive  approaches.  In  our  new  schemes  the 
reduction  is  achieved  by  using  a  priori  information  that  is  not  yet  being  exploited.  In  this 
report,  we  focus  on  using  our  knowledge  on  the  structure  of  ground  clutter,  but  other  prior 
knowledge  could  also  be  exploited.  For  example,  one  might  want  to  incorporate 
knowledge  of  jammers  or  interference  from  previous  scans  or  one  might  have  knowledge 
of  previously  detected  targets.  Using  prior  knowledge  should  allow  our  new  schemes  to 
perform  well  even  when  given  only  small  amounts  of  reference  data,  which  may  have  a 
statistical  description  which  is  slightly  different  from  that  for  the  cell-under-test. 

A.  Traditional  schemes 

Assume  that  the  observations  to  be  processed  are  taken  from  M  different  pulse  returns 

received  by  N  antenna  elements  from  the  k  range  cell.  Each  return  is  assumed  to 
contain  a  possible  signal  in  additive  noise-plus-clutter.  Denote  the  observation 
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corresponding  to  the  i"'  pulse  at  the  antenna  element  as  X.j.  Each  observation  is  a 

complex  number  corresponding  to  the  in-phase  and  quadrature  components  of  the  received 
waveform.  The  observations  to  be  processed  are  ordered  as 

X 1^  »  ^21  ’  * "  *  ’ ’ '^12  ’  * "  ’ 

where  X  denotes  the  transpose  of  the  vector  X . 

A  reasonably  large  set  of  STAP  processing  schemes  can  be  described  as  follows.  The 
description  is  broken  into  three  parts;  pre-processing,  adaptive  processing  and  post¬ 
processing.  The  pre-Processing  can  be  described  by 

X,(p)  =  {A^®BXX,.  p  =  0,1,2 . P-l  (2) 

where  is  given  in  (1).  In  (2)  and  are  scheme-dependent  matrices.  The 
adaptive  processing  can  be  described  by 

y,(p)  =  S''/?;‘(P)X.(p)/1>,.  p  =  0,1,2,...,  p-l  (3) 

where  S  =  ®  (4) 

is  the  scheme-dependent  temporal  steering  vector  and  is  the  scheme-dependent 

spatial  steering  vector.  In  (4)  (S)  stands  for  an  outer  product.  The  quantity  is  the 

normalization  needed  to  provide  a  constant  false  alarm  rate  (CFAR)  in  homogeneous 
clutter  and  is  given  by 


The  post-processing  is  not  necessary  if  (p)  is  the  final  output  of  interest.  In  this  case, 
\yk  (P)l  compared  to  a  threshold  directly  to  decide  if  target  is  present.  If  yi^ip)  is  not 
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the  final  output  of  interest,  we  should  assemble  all  the  outputs  into  a  vector  as 
^  .then  compute 

(6) 

typically,  is  the  m'*  column  of  a  P  X  P  filter  matrix  P ,  and  Z*  „  is  the  final  output 
whose  magnitude  is  compared  to  a  threshold  to  make  a  decision. 

In  the  traditional  schemes  we  considered,  the  estimated  covariance  matrix  is  calculated  as 

1  k+LI2+\ 

R,ip)=-  (7) 

by  averaging  over  L  range  cells  (assuming  L  is  even)  surrounding  the  k'’'  range  cell 
excluding  the  cell-under-test  and  the  two  closest  range  cells. 

B.  Modified  schemes 


In  order  to  use  our  prior  knowledge  of  the  nature  of  ground  clutter,  namely  the  known 
structure  of  the  clutter  ridge  in  angle-Doppler  space,  we  use  a  simple  model  for  ground 
clutter  proposed  in  [1].  Assume  that  the  clutter  portion  of  the  reference  samples  is  Gaussian 
distributed  with  the  two-dimensional  power  spectrum  density  (psd)  described  in  [1] 

f  i  f*  n  /  /•  y,  \2  ^ 


d 


c4 


exp 


{f.-fcJ)  .{f.-fca) 


2a 


ft  4 


(8) 


which  is  a  function  of  normalized  Doppler  frequency  and  spatial  frequency  .  The 
psd  in  (8)  consists  of  L  Gaussian-shaped  humps,  the  d  th  of  which  is  centered  at 
(/f  ^fs)~  {fct^  ’  fcs^  )  amplitude  controlled  by  and  a  spread  in  angle  and 
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2  2 

Doppler  controlled  by  (7^^  ^  and  (J^  ^ .  Our  approach  will  assume  the  clutter  samples  are 

described  by  the  psd  in  (8)  and  the  data  will  be  used  to  estimate  some  of  the  parameters  in 
(8).  The  other  parameters  in  (8)  will  be  fixed  at  values  chosen  to  fit  a  wide  set  of  measured 
data.  Once  the  estimates  for  the  parameters  in  (8)  are  found,  the  covariance  matrix  for 
clutter  described  by  (8)  is  given  by  [1] 

d=l 

where  ^  and  ^  are  Toeplitz  matrices  specified  by 

-2{nafi  d)^+i2;r«.d 

6  '"C 

and 

_  fr^  -2{itak4)^+i2Tfcs,d  -2.(7cafs4(N-\))^+i(N-\)2T^cs,d 

=Toeplitz<  1  e 
respectively. 

Now  we  can  modify  the  traditional  schemes  by  using  the  covariance  matrix  in  (9), 
possibly  with  slight  modification,  to  replace  the  covariance  matrix  estimate  of  (7). 

For  our  modified  schemes,  the  covariance  matrix  is  derived  as 

R,  ip)  =  e{(x,  (p)  -  E{X.  (p)})(x,  (p)  -  e{x.  (p)})"} 

=  (A,  ®  ®  B,) 

d=\ 


=  Toeplitz  I  1 
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(12) 


d^\ 

where  (9)  was  used  along  with  Theorem  8.8.2  and  Theorem  8.8.6  in  [61.  C  .  and  C  . 
are  given  in  (10)  and  (11). 

To  illustrate  the  modified  approach  just  described,  the  Sample  Matrix  Inversion  algorithm 
(SMI),  the  Adaptive  Displaced  Phase-Centered  Antenna  algorithm  (ADPCA),  the 
Extended  Factored  Post-Doppler  algorithm  (EFA)  and  the  Joint-Domain  Localized 
algorithm  (JDL)  are  used  as  examples.  We  call  the  new  schemes  modified  SMI,  modified 
ADPCA,  modified  EFA  and  modified  JDL  respectively. 

For  the  modified  schemes,  we  do  not  develop  our  covariance  matrix  estimate  by 
performing  averages  over  the  adjacent  range  cells.  Instead  the  covariance  matrix  is 
determined  by  the  model  parameters  in  (9)  which  are  obtained  from  the  range  cell  under 
test.  Thus  mismatch  between  the  statistics  of  the  reference  data  and  the  data  from  the  cell 
under  test  is  avoided,  which  can  be  a  problem  for  traditional  schemes  in  a  highly  non- 
homogeneous  environment.  In  the  traditional  SMI,  ADPCA,  EFA  and  JDL  schemes  used 
in  our  comparison,  we  use  a  set  of  reference  samples  on  either  side  of  the  cell-under-test. 
In  our  examples,  66  cells  (33  on  either  side  of  the  cell-under-test)  are  used  for  SMI, 
^‘^PCA  and  EFA.  10  cells  are  used  for  JDL  (5  on  each  side  of  the  cell-under-test).  A 
covariance  matrix  is  formed  by  averaging  the  covariance  matrix  estimated  from  each  of 
those  adjacent  cells. 

Obviously,  the  parameters  for  the  model  in  (8)  determine  the  required  covariance  matrix. 
Our  tests  indicate  we  do  not  need  to  estimate  all  these  parameters  for  every  range  cell  to 
achieve  good  performance.  Our  approach  assumes  some  of  the  parameters  are  fixed  at 
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reasonable  values  previously  estimated.  In  particular,  the  parameters  fixed  are  observed  to 
have  little  impact  on  performance.  We  only  devote  computational  resources  to  obtain  the 
other  ‘key’  parameters.  Below  is  the  discussion  of  each  of  the  ‘key’  parameters  and  how  to 
choose  proper  values  for  them. 


3.  Parameter  choice 


We  use  data  that  comes  from  the  MCARM  [7,8]  database  flight  5,  acquisition  575  in  our 
discussion.  In  order  to  obtain  an  estimate  of  the  psd  in  the  angle-Doppler  domain,  we 
apply  an  FFT  to  transform  to  the  Doppler  domain  first,  then  we  use  the  normalized  steering 
vectors  provided  with  the  MCARM  database  to  transform  to  the  angle-Doppler  domain. 
This  would  be  essentially  equivalent  to  performing  a  two  dimensional  FFT  if  the  antenna 
array  were  a  uniformly  spaced  linear  array. 

A.  Slope  of  clutter  ridge 

The  most  important  parameter  of  the  model  in  (8)  is  the  slope  of  the  clutter  ridge  along 
which  we  center  the  Gaussian  functions.  The  slope  of  clutter  ridge  is  related  to  airplane’s 
velocity  which  might  be  known  directly  from  instruments  on  the  aircraft.  In  our 
experiments,  we  have  used  some  fast  image  processing  algorithms  to  find  this  slope 
directly  from  data. 
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B.  Number,  amplitudes,  and  positions  of  humps 


For  simplicity,  we  fix  the  number  of  humps  used  and  place  them  uniformly  along  the 
clutter  ridge.  If  the  number  of  humps  is  too  small,  the  model  cannot  represent  the  real 
clutter  well.  On  the  other  hand  if  the  number  of  humps  is  too  large,  computational 
complexity  increases  with  little  performance  improvement.  Extensive  testing  indicated 
that  using  4  or  5  humps  is  a  good  choice. 

We  also  fix  the  amplitude  of  the  humps.  We  found  that  reasonably  good  performance  is 
obtained  by  setting  the  amplitude  of  each  hump  to  be  approximately  15%  greater  than  the 
maximum  value  of  the  estimated  psd.  In  this  report,  all  the  results  shown  use  this  choice. 

C.  :  parameters  controlling  the  spread  of  the  hump 

Extensive  experiments  showed  that  spread  parameters  should  be  chosen  differently  for  the 
pre-Doppler  schemes  (such  as  the  modified  SMI  and  the  modified  ADPCA)  and  the  post- 
Doppler  schemes  (such  as  the  modified  EFA  and  the  modified  JDL). 

For  the  pre-Doppler  schemes,  two  ways  are  suggested  to  set  the  parameters  <7^  ^ ^  • 

One  is  setting  a  ‘fixed’  value  for  all  humps.  The  fixed  value  should  be  set  larger  than  that 
the  estimated  psd  implies.  The  other  method  attempts  to  provide  values  that  more  closely 
match  the  estimated  psd.  In  order  to  explain  when  each  of  these  methods  should  be  used 
we  show  some  examples.  We  inserted  a  single  synthetic  moving  target  with  a  specified 
normalized  Doppler  frequency  into  a  specified  range  cell  to  compare  the  detection 
performance  for  the  two  cases. 
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We  present  results  which  show  the  magnitude  of  the  test  statistic  for  the  range  cell-under¬ 
test  and  for  surrounding  range  cells  as  in  [5].  If  the  response  at  the  range  cell  where  the 
target  is  inserted  is  much  larger  than  surrounding  range  cells,  then  we  judge  performance  as 
being  very  good.  We  normalized  all  adaptive  weights  by  a  scalar  (see  (3)  in  [5]),  which 
provides  a  CFAR  test  for  homogeneous  cases.  The  amplitude  of  the  target  is  set  at  0.05. 

In  the  first  example,  we  set  (7^^.  =  (7^^.  =  0.1,  for  i  —  and  we  inserted  a 

single  target  with  spatial  frequency’  0.164  and  normalized  Doppler  frequency  0.078,  0.156, 
or  0.312  in  range  bin  350.  In  a  second  set  of  results  we  inserted  a  single  target  in  range  bin 
415  at  the  same  set  of  spatial  and  Doppler  frequencies.  We  present  the  results  for  modified 
ADPCA.  The  results  are  illustrated  in  Fig.  1. 

In  the  second  example,  where  we  choose  the  spreads  to  more  closely  match  what  is 
observed  in  the  psd  estimate,  we  set  C7^ .  =  .  =  0.01,  i  =  1,5,  and 

(7^.  =0.02,  /  =  2,3,4.  For  this  case,  we  get  the  performance  shown  in 

Fig.  2. 

We  note  that  for  the  target  with  normalized  Doppler  frequency  0.078  (whatever  the  range 
bin),  the  performance  in  the  second  case  is  much  better  than  the  first  one.  On  the  other 
hand,  for  the  target  with  normalized  Doppler  frequency  0.156  and  0.312  (whatever  the 

range  bin) ,  the  first  case  with  ‘fixed’  and  ‘larger’  C7^  ^ ^  will  be  better.  The  reason  is 

that  when  target  has  the  normalized  Doppler  frequency  0.078,  its  position  is  near  the  clutter 
ridge.  If  the  spread  parameters  are  too  large,  the  target  itself  will  be  significantly 

'  In  all  the  tests  of  this  report,  the  spatial  frequency  of  target  is  fixed  at  0.164.  Experiments  showed  that  these 
results  are  representative. 
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suppressed.  Thus,  in  this  case,  we  should  choose  the  parameters  ‘carefully’  and  use 
parameters  which  match  the  true  parameters  well.  In  the  other  cases,  when  the  target  has 
the  normalized  Doppler  frequency  0.156  or  0.312,  the  target  position  is  not  so  close  to  the 
clutter  ridge.  In  this  case  it  is  better  to  use  larger  spread  parameters  to  provide  increased 
clutter  suppression,  since  the  target  will  not  be  suppressed  when  this  is  done.  In  practice, 
we  can  compare  the  relative  target  to  clutter  ridge  position  to  a  threshold  to  decide  which 

way  the  parameters  should  be  set.  While  the  examples  provided  here  are  for 

modified  ADPCA,  further  tests  demonstrated  that  the  conclusion  is  the  same  for  all  the 
other  pre-Doppler  processing  schemes  we  tested. 

For  the  post-Doppler  schemes,  we  found  that  no  matter  where  the  target  is  in  the  angle- 
Doppler  domain,  the  best  choice  for  the  spread  parameters  (7^  ^ ,  (7^^  ^  is  using  those  that 
closely  match  the  observed  psd. 

Four  modified  schemes  are  analyzed  in  next  section  and  parameters  are  chosen  according 
to  the  discussion  in  this  section. 


4.  Performance  comparison 

Here,  we  will  investigate  the  performance  obtained  when  using  the  modified  SMI,  the 
modified  ADPCA,  the  modified  EFA  and  the  modified  JDL  algorithms.  Specifically  we 
compare  the  performance  of  the  modified  scheme  with  that  for  the  corresponding 
traditional  scheme.  The  comparisons  are  similar  to  those  in  Fig.  1  and  Fig.  2. 
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A.  SMI 


In  our  tests  of  SMI  we  use  M=22  pulses  and  N=22  antennas  in  the  adaptive  processing. 
For  the  SMI  scheme,  we  have  =  /^,  in  (2),  where  are  MXM 

and  N  X  N  identity  matrices  respectively,  and  are  the  target  temporal  and  spatial 
response  vectors  respectively  [1].  Post-processing  is  not  necessary  for  SMI. 

First,  we  inserted  a  single  target  with  spatial  frequency  0.164  and  with  normalized 
Doppler  frequency  of  either  0.078,  0.156,  or  0.312  in  range  bin  350.  In  a  second  set  of 
tests,  a  single  target  with  spatial  frequency  0.164  and  with  normalized  Doppler  frequency 
of  either  0.078,  0.156,  or  0.312  was  inserted  in  range  bin  415.  In  each  case  the  amplitude 
at  the  target  was  0.05.  Using  modified  SMI,  we  obtain  the  results  in  Fig.  3.  Traditional 
SMI  yielded  the  results  in  Fig.  4.  For  the  case  of  range  bin  350,  the  performance  of  the  two 
schemes  is  close  except  when  the  target  is  far  from  the  clutter  ridge.  In  particular,  the 
performance  of  the  traditional  SMI  scheme  in  Fig.  3(e)  is  much  better  than  that  of  the 
modified  scheme  in  Fig.  4(e).  On  the  other  hand,  the  modified  scheme  is  always  much 
better  for  range  bin  415  as  illustrated  in  Fig.  3  and  Fig.  4. 

In  order  to  explain  the  difference  between  the  results  for  range  bin  350  and  those  for  range 
bin  415,  we  consider  Fig.  5  which  shows  average  energy  as  a  function  of  range  bin.  Near 
range  bin  350,  the  variation  is  not  too  great.  In  this  case  it  is  advantageous  to  average  over 
range  when  estimating  the  covariance  matrix  and  this  leads  to  the  good  performance  of 
traditional  SMI.  Since  traditional  SMI  is  optimum  when  the  estimate  of  the  covariance 
matrix  is  perfect,  it  is  reasonable  that  we  cannot  outperform  SMI  in  some  cases  where  this 
is  close  to  being  true.  On  the  other  hand.  Fig.  5  shows  a  large  jump  in  energy  near  range 
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bin  400.  So  this  would  impact  the  averaging  performed  in  the  traditional  schemes  when  a 
target  is  inserted  in  range  bin  415.  This  non-homogeneous  reference  data  leads  to  the  poor 
performance  of  traditional  SMI  and  in  turn  to  the  superior  performance  of  modified  SMI. 


B.  ADPCA 


For  ADPCA,  the  matrices  A  and  B  used  in  (2)  are 

p  p 


P  k, 


p  N 


where  the  notation  0^^^  refers  to  a  qX  m  matrix  of  zeros.  In  this  research  M  is  set  to  128, 
N  is  set  to  22  and  is  set  to  3.  The  A^  in  (13)  groups  the  M  pulses  in  subgroups  of  size 
3  with  each  group  overlapping  by  2  pulses.  In  the  adaptive  processing  part, 
•S',  =  [l  —  2  l] ,  and  is  the  target  spatial  response  vector.  In  the  Post-processing 

fm  is  the  Doppler  filter  corresponding  to  the  target  Doppler  frequency. 


There  are  a  number  of  advantages  to  partially  adaptive  schemes  in  general  and  ADPCA  in 
particular,  see  [3,  4]  for  a  discussion  on  these  issues.  In  order  to  compare  the  performance 
of  the  traditional  ADPCA  scheme  and  the  modified  ADPCA  scheme,  we  consider  the  same 
cases  as  for  SMI.  In  addition,  we  also  considered  targets  with  normalized  Doppler 
frequency  of  0.039,  0.117  or  0.234.  The  results  for  modified  ADPCA  and  traditional 
ADPCA  are  shown  in  Fig.  6  and  Fig.  7  respectively.  The  modified  scheme  is  clearly  better 
for  cases  where  the  target  has  the  normalized  Doppler  frequency  0.039  or  0.078,  (when  the 
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target  is  near  the  clutter  ridge).  For  cases  where  the  target  has  the  normalized  Doppler 
frequency  0.1 17  or  0.156,  (where  the  clutter  ridge  is  not  so  close  to  the  target),  the  modified 
scheme  results  are  significantly  better  for  range  bin  415.  However,  performance  is  about 
the  same  when  the  target  is  in  range  bin  350.  For  cases  where  the  target  has  the  normalized 
Doppler  frequency  0.234  or  0.312,  (where  the  clutter  ridge  is  far  from  the  target),  the 
traditional  scheme  is  better  for  range  bin  350  while  the  schemes  provide  about  the  same 
performance  for  range  bin  415.  Taken  collectively  these  results  indicate  the  modified 
scheme  should  be  used  in  non-stationary  cases.  The  advantages  seem  to  be  especially  clear 
when  the  target  is  near  the  clutter  ridge. 


C.  EFA 

For  the  EFA  scheme,  the  pre-processing  is  as  in  (2)  with 

A  =  [4-.  -  4 . 4«  ]  =  I  (14) 

where  f  ,  f  are  three  adjacent  Doppler  filters  and  corresponds  to  the 

normalized  Doppler  frequency  of  target.  The  adaptive  processing  uses  =  [0  1  0]^ 

and  set  to  the  target  spatial  response  vector.  Only  processing  for  a  single  p  must  be 
performed  and  so  post-processing  is  not  needed. 

The  results  for  modified  EFA  and  traditional  EFA  are  given  in  Fig.  8  and  Fig.  9.  The 
results  for  range  bin  350  indicate  that  when  the  target  is  inserted  far  from  the  clutter  ridge 
(target  with  normalized  Doppler  frequency  0.234  or  larger),  the  modified  scheme  achieves 
significantly  better  performance.  In  the  other  cases,  the  modified  scheme  is  slightly  better. 
For  range  bin  415,  the  modified  scheme  provides  better  results  in  all  cases  except  the  case 
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in  which  the  target  has  normalized  Doppler  frequency  0.117.  Since  we  have  already 
identified  range  bin  350  as  an  almost  stationary  case  and  range  bin  415  as  a  highly  non¬ 
stationary  case,  these  results  are  reasonable.  It  is  known  that  in  a  nearly  homogeneous 
case,  we  should  use  the  traditional  scheme  since  in  this  case  the  estimate  of  the  covariance 
matrix  is  a  maximum  likelihood  estimate.  In  a  highly  non-homogeneous  case,  it  is 
reasonable  that  the  modified  scheme  can  provide  better  performance. 


D.  JDL 

For  JDL  the  matrices  in  (2)  should  be  taken  as 

~  ’ /p  ’ -/p+l  ]  ’  ^  p  Sqi  Sq+\\^  (15) 

where  /^_j ,  are  three  Doppler  filters  and  is  that  filter  which  corresponds  to 

the  normalized  Doppler  frequency  of  target.  j  and  are  Doppler  filters 

corresponding  to  neighboring  Doppler  frequencies.  are  three  spatial  filters 

and  g^  is  that  filter  which  corresponds  the  spatial  steering  vector  of  the  target.  g^_^  and 
Sq+\  spatial  filters  corresponding  to  neighboring  spatial  frequencies.  The  steering 
vector  S  in  the  adaptive  processing  is 

S  =  S,(8)5,=[0  1  0]®[0  1  0]  (16) 

As  for  EFA,  processing  must  be  performed  for  only  one  p  so  post-processing  is  not 
needed. 
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We  repeat  the  tests  made  for  EFA  for  JDL  to  compare  the  performance  of  the  modified 
and  the  traditional  schemes.  The  results  are  given  in  Fig.  10  and  Fig.  11  respectively. 
When  the  target  is  inserted  in  range  bin  350,  the  performance  of  modified  scheme  is  much 
worse  than  the  traditional  scheme.  The  problem  seems  to  be  that  traditional  JDL  performs 
very  well  in  these  cases.  When  the  target  is  inserted  in  range  bin  415  far  from  the  clutter 
ridge  (target  with  normalized  Doppler  frequency  0.234  or  0.312),  the  modified  scheme  is 
better.  In  the  other  cases  considered  for  range  bin  415,  the  performance  is  about  the  same. 
Thus  the  performance  of  modified  JDL  is  rather  poor  in  the  homogeneous  case.  However 
in  most  non-homogeneous  cases,  it  performs  well. 


5.  Conclusion 


We  have  suggested  the  use  of  some  modified  schemes  for  adaptive  airborne  radar,  which 
use  the  knowledge  of  the  structure  of  the  interference.  The  specific  numerical  results  given 
use  a  model  which  only  assumes  ground  clutter,  while  these  results  are  obtained  for  a 
measured  radar  data  set  which  includes  other  interference. 

From  the  comparisons,  it  appears  that  the  modified  scheme  can  provide  good  performance 
in  the  cases  that  are  traditionally  considered  to  be  difficult.  It  performs  better  in  non- 
homogeneous  cases  when  compared  to  schemes  which  average  over  range  to  estimate  the 
required  covariance  matrix.  This  is  reasonable  since  the  modified  schemes  obtain  the 
covariance  matrix  from  a  single  range  cell.  However,  when  the  target  is  inserted  in  a 
homogeneous  environment,  in  most  cases  the  traditional  schemes  are  better. 
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e)  range  bin  350,  Doppler  frequency  0.3 12  f)  range  bin  415,  Doppler  frequency  0.3 12 
Figure  1)  Performance  comparison  for  target  inserted  in  different  range  bins  and  with 
several  different  normalized  Doppler  frequencies  for  modified  ADPCA  when 


^ ^ fs4  1,...^. 
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a)  range  bin  350,  Doppler  frequency  0.078  b)  range  bin  415,  Doppler  frequency  0.078 


e)  range  bin  350,  Doppler  frequency  0.3 12  f)  range  bin  415,  Doppler  frequency  0.3 12 
Figure  2)  Performance  comparison  for  target  inserted  in  different  range  bins  and 
with  several  different  normalized  Doppler  frequencies  for  modified  ADPCA  when 

=  0.01 ,  i  =  u,  and  cr  =  C7  =  0.02 ,  ;•  =  2,3,4 . 


330  340  3S0  380  370  380 


a)  range  bin  350,  Doppler  frequency  0.078  b)  range  bin  415,  Doppler  frequency  0.078 


A  M  lj'< 


c)  range  bin  350,  Doppler  frequency  0. 156  d)  range  bin  415,  Doppler  frequency  0. 156 


e)  range  bin  350,  Doppler  frequency  0.312  f)  range  bin  415,  Doppler  frequency  0.312 


Figure  3)  Performance  for  modified  SMI  when  target  is  inserted  in  a  particular  range 
bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized  spatial 
frequency  of  the  target  is  0.164. 


e)  range  bin  350,  Doppler  frequency  0.312  f)  range  bin  415,  Doppler  frequency  0.312 


Figure  4)  Performance  for  traditional  SMI  when  target  is  inserted  in  a  particular  range 
bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized  spatial 
frequency  of  the  target  is  0. 164. 
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Figure  5)  Average  energy  as  a  function  of  range  bin 
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k)  range  bin  415,  Doppler  frequency  0.234  1)  range  bin  415,  Doppler  frequency  0.312 
Figure  6)  Performance  for  modified  ADPCA  when  target  is  inserted  in  a  particular 
range  bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized 
spatial  frequency  of  the  target  is  0.164. 
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e)  range  bin  350,  Doppler  frequency  0.234  f)  rangebin350,  Doppler  frequency  0.31! 


Figure  7)  Performance  for  traditional  ADPCA  when  target  is  inserted  in  a  particular 
range  bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized 
spatial  frequency  of  the  target  is  0.164. 
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390  400  410  420  430  440  450  460  390  400  410  420  430  440  490  460 

k)  range  bin  415,  doppler  frequency  0.234  1)  range  bin  415,  doppler  frequency  0.312 

Figure  8)  Performance  for  modified  EFA  when  target  is  inserted  in  a  particular  range 
bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized  spatial 
frequency  of  the  target  is  0. 164. 
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g)  range  bin  415,  doppler  frequency  0.039  h)  range  bin  415,  doppler  frequency  0.078 


i)  range  bin  415,  doppler  frequency  0.117  j)  range  bin  415,  doppler  frequency  0.156 


k)  range  bin  415,  doppler  frequency  0.234  1)  range  bin  415,  doppler  frequency  0.312 

Figure  9)  Performance  for  traditional  EFA  when  target  is  inserted  in  a  particular 
range  bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized 
spatial  frequency  of  the  target  is  0.164. 
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i)  range  bin  415,  doppler  frequency  0. 1 17  j)  range  bin  415,  doppler  frequency  0. 156 


Figure  10)  Performance  for  modified  JDL  when  target  is  inserted  in  a  particular  range 
bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized  spatial 
frequency  of  the  target  is  0.164. 


k)  range  bin  415,  doppler  frequency  0.234  1)  range  bin  415,  doppler  frequency  0.312 


Figure  11)  Performance  for  traditional  JDL  when  target  is  inserted  in  a  particular 
range  bin  with  a  particular  normalized  Doppler  frequency.  In  each  case  the  normalized 
spatial  frequency  of  the  target  is  0.164. 


71 


Chapter  3  :  Performance  Analysis 
Applicable  to  both  Tasks  4.1.1  and  4.1.2 


Analysis  of  STAR  Algorithms  for  Cases  with  Mismatched  Statistics 
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Introduction 


In  most  adaptive  radar  implementations,  the  clutter-plus-noise  in  the  cell-under-test  is 
characterized  using  samples  taken  from  range  cells  that  neighbor  the  cell-under-test.  This 
can  lead  to  a  mismatch  between  the  true  clutter-plus-noise  statistics  (in  the  cell-under¬ 
test)  and  those  used  to  design  the  adaptive  processing  scheme.  Such  mismatches  can 
occur  in  nonhomogeneous  noise-plus-clutter  cases  [1].  The  purpose  of  this  research  is  to 
develop  analytical  formulas  that  characterize  the  loss  in  performance  due  to  this 
mismatch.  These  formulas  should  be  useful  as  a  performance  metric  for  developing 
STAP  schemes  that  are  not  so  sensitive  to  mismatch.  The  focus  is  on  the  constant  false 
alarm  rate  (CFAR)  version  of  the  sample  matrix  inversion  (SMI)  algorithm,  but  the 
analysis  is  extended  to  other  algorithms  also.  The  noise-plus-clutter,  which  includes 
clutter,  jamming,  and  noise,  is  assumed  to  be  a  complex  Gaussian  process. 

Consider  a  CFAR  version  of  the  SMI  algorithm  of  the  type  discussed  in  [2,3]  with  the 
test  statistic 

(1) 

In  (1),  the  observed  N  dimensional  complex  vector  x  consists  of  either  zero-mean 
complex  Gaussian  noise-plus-clutter  with  covariance  matrix  /?,  or  signal  plus  zero-mean 
complex  Gaussian  noise-plus-clutter  with  covariance  matrix  R,.  The  signal  which  is 
added  to  the  noise-plus-clutter  is  k  s  where  5  is  a  unit  length  signal  vector  which  is 
completely  known  and  /cis  an  unknown  complex  constant  The  magnitude  of  vsets  the 
signal-to-noise  ratio.  The  denominator  of  (1)  provides  the  correct  normalization  for 
CFAR  (for  cases  without  mismatch).  While  it  would  be  desirable  to  use  the  true 
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covariance*  R,  in  (1)  in  place  of  Re,  Ri  is  not  available.  Instead  an  estimated  covariance 
matrix  Rg  is  used.  Re  is  obtained  from  the  maximum  likelihood  (ML)  estimate  for  the  case 
where  a  set  of  independent  reference  vectors  x(k),  k  =  are  available,  all  with  the 

same  distribution  as  the  cell-under-text  data  x.  Specifically,  Rg  is  taken  as 

R,  ='£,x(k)x{k)"  (2) 

ifc=l 

which  differs  from  the  ML  estimate  by  a  scale  factor.  Any  scale  factor  modification  of 
(2),  and  thus  Rg,  leads  to  the  inverse  of  this  same  scale  factor  change  in  (1)  which  can  be 
incorporated  into  the  threshold  (1)  is  compared  to.  Note  that  in  practice  one  frequently 
finds  Rg  is  not  close  to  a  scaled  version  of  Rt,  due  to  a  variety  of  reasons.  In  fact,  the 
expected  value  of  Rg  may  differ  from  L  Rt,  which  implies  that  even  for  very  large  L,  the 
value  of  (l/L)Rg  will  not  converge  to  Rt.  For  the  purpose  of  this  paper  we  assume  a 
mismatch  in  the  reference  data  such  that  (l/L)E{Rg}  -  Rsd^Rt-  For  simplicity  we  assume 
that  the  reference  data  vectors  are  independent  and  identically  distributed  (iid)  with  a 
zero-mean  complex  Gaussian  distribution  with  covariance  matrix  R^d.  Also,  the 
reference  data  are  independent  from  the  data  from  the  cell-under-test. 


Distribution  of  the  Test  Statistic 

First  apply  a  coordinate  transform,  which  consists  of  multiplication  by  Rsd'^^,  to  the 
observed  vector  from  the  cell-under  test,  the  reference  data  vectors  and  to  the  signal 
vector.  This  transform  whitens  the  reference  data.  Next,  normalize  the  transformed 
signal  vector  so  that  it  is  again  a  unit  vector.  Then,  call  this  transformed  unit  signal 
vector  V.  The  operation  is  equivalent  to  scaling  the  signal  Rsd^^s  by  dividing  by 
(s^Rsd  So  that  signal-to-noise  ratio  is  maintained,  the  transformed  signal  is  taken  to 
he  fiv  =  [(s^Rsd  k]  V.  The  important  parameter  Ij3l^  is  called  the  signal-to-secondary 

‘  This  would  lead  to  (1)  being  optimum  [9], 
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noise  ratio  (SSNR).  It  plays  a  role  similar  to  the  role  the  signal-to-noise  ratio  plays  in 
cases  without  mismatch.  Next,  another  transformation  of  coordinates,  which  consists  of 
multiplication  by  (  v,  f,  is  made  to  the  cell-under-test  vector,  the  reference  data 
vectors  and  to  the  signal  vector.  Here  By  is  a  matrix  whose  N-1  rows  consist  of  a  set  of  N 
dimensional  vectors  which  span  the  space  orthogonal  to  v.  The  overall  matrix  (  v,  B/^  f 
is  taken  to  be  unitary  so  that  the  signal  and  noise  powers  are  preserved.  Call  the 
transformed  vector  for  the  cell-under-test  y  =  (  d,  b^f.  Here  is  a  scalar  that  describes 
the  component  of  the  vector  Rsd  ^^x  which  lies  in  the  direction  of  Rsd  ’^s.  On  the  other 
hand,  b  is  an  N-1  dimensional  vector  describing  the  component  of  Rsd^x  that  is 
orthogonal  to  Rsd'^^s.  For  the  transformed  reference  vectors  we  employ  the  notation  y(k) 
=  (  d(k),  b(kff.  The  transformations  made  have  a  similar  flavor  to  those  used  in  [4]  for 
a  similar  purpose.  These  transformations  also  convert  the  adaptive  filtering  structure 
from  the  standard  “direct-form”  to  the  “generalized  sidelobe  canceler”  form  [5,6].  For 
generality,  we  could  assume  that  the  b  and  b(k)  components  are  multiplied  by  the  matrix 
m"  to  obtain  z  =  u^b  and  z(k)=  u%(k)  and  that  z  and  z(k),  k-l,...,L  are  used  to  form  the 
test  statistic.  This  allows  consideration  of  rank  reduction  schemes.  For  cases  without 
rank  reduction,  just  take  m  to  be  an  identity  matrix.  If  rank  reduction  is  desired  then  u  is 
taken  to  be  a  column  rank  M  unitary  matrix  where  M  <  N-1.  Under  the  assumptions 
outlined,  we  find  that  (1)  becomes  (see  Appendix) 


where 


(3) 


f^=iz{k)d{kf  (4) 

*=i 

d]  =  '^d(k)dik)''  (5) 

k=l 
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and 


R^='^z(^k)zik)''  (6) 

k=l 

Using  the  transformations  described  yields 

E{RJ  =  R^  =  =  Lu^u 

Note  that  without  rank  reduction  (u  is  an  identity)  —  LI,  where  I  is  an  identity  matrix. 
Now  condition  on  z,  z(l),...,z(L)  and  note  that  the  quantity  inside  the  absolute  value  signs 
in  the  denominator  of  (3)  can  be  expressed  in  matrix  form  as  (as  in  [6]  for  the  case 
without  mismatch) 

DID"  =  d{i  -  Z"R-'z)d"  (7) 


where 


D  =  id(,l),d(2),...,d{L)) 

and 

Z  =  (z(l),z(2),...,z(L)) 

It  is  possible  to  demonstrate  (see  Appendix)  that  the  matrix  Z  in  (7)  has  rank  L-M  where 
M^-1  (M=N-l  in  cases  with  no  rank  reduction).  Recall  M  is  the  column  rank  of  the 

rank  reduction  matrix  u.  It  is  also  possible  to  demonstrate  (see  Appendix)  that  R~^Z 
is  an  idempotent  projection  matrix.  Using  this  and  results  from  [7,  pp.  113-118]  and  [8, 
pp.  27-29],  it  follows  that  the  denominator  of  (3),  when  conditioned  on  z,  z(l),...z(L),  is  a 
constant  factor  of  V^  times  a  central  chi-square  distributed  random  variable,  with  2{L-M) 
degrees  of  freedom.  Such  a  random  variable  is  denoted  by  In  a  similar  way  the 

term  in  the  numerator  of  (3),  inside  the  I  P  is 
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(8) 


y,^d-J^d{k)z{k)^R:h 

k=\ 

which  is  zero-mean,  complex  Gaussian  distributed  when  conditioned  on  z,  z(l),..-z(L)  and 
assuming  signal  is  not  present.  Under  the  same  conditioning,  when  signal  is  present,  yr  is 
complex  Gaussian  distributed  with  mean  p.  Under  either  signal  absent  or  signal  present 
the  variance  of  yr  (same  conditioning)  is 

Var{y^}=  —  =  Var{d}-{-  z^R'^z 

P 

(see  Appendix  for  justification  and  the  value  of  Var{d}).  Next,  we  properly  normalize 
(multiply  by  p)  the  numerator  of  (3)  so  it  is  the  square  of  a  unit-variance  complex 
Gaussian.  Conditioned  on  z,  z(l),...z(L),  another  Theorem  from  [7,  pp.  113-118]  shows 
that  1  yr\^  is  a  constant  factor  of  V2  times  a  non-central  chi-square  distributed  random 
variable  with  2  degrees  of  freedom,  X2^(p  when  signal  is  present.  When  signal  is 
absent,  the  same  holds  true  with  P=0.  Thus,  in  sununary,  conditioned  on  z,  z(l),.-.z(L) 
the  test  statistic  in  (3)  is  the  ratio  of  a  non-central  chi-squared  random  variable  with  2 
degrees  of  freedom  to  p  times  a  central  chi-squared  random  variable  with  2(L-M) 

degrees  of  freedom 


Under  the  assumption  that  the  reference  data  are  independent  from  the  data  from  the 

cell-under-test,  the  two  chi-squared  random  variables  in  the  numerator  and  denominator 
of  Aamf  are  independent.  Then,  from  [8,  p.  24-25]  (with  slight  extension  for  complex 
quantities)  the  probability  density  function  (pdf )  of  the  test  statistic  is 
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/ (^Ip)  - 


(L-M)pexp(-p|j3'l) 


(l  +  pA) 


L-M+1 


.F,(L-M  +  1;1;(2) 


(9) 


where 


\gyx 

1 4*  pA 


Equation  (9)  gives  the  pdf  of  Aamf  when  conditioned  on  z,  z(l),...z(L)  which  is 
equivalent  to  conditioning  on  p  as  indicated  by  the  notation  used  in  (9).  Further  the 
result  in  (9)  was  obtained  under  the  signal  present  (Hj)  hypothesis.  In  (9)  iFi(a;b;c)  is 
used  to  denote  the  generalized  hypergeometric  function  [8,  p.  20].  The  result  in  (9)  is 
valid  for  signal  absent  (Ho)  if  one  sets  P=0. 

To  find  the  unconditional  pdf  of  Aamf,  (9)  must  be  averaged  over  the  random  fluctuations 
due  to  estimating  the  covariance  matrix  in  (3).  It  is  interesting  that  the  loss  due  to  this 
estimation  enters  (9)  only  through  p.  It  is  helpful  to  notice  that  p  is  directly  related  to  the 
observations,  as  stated  previously,  by 


P  = 


l  +  z^'R^'z 


(10) 


where 


Z 


Thus  the  pdf  for  p  has  support  for  p  between  0  and  1/Var{d}  (otherwise  the  pdf  is  zero). 
To  perform  the  required  averaging,  first  note  that  the  estimated  covariance  matrix  in  (2) 
is  characterized  by  the  Wishart  distribution  [8,9]  when  the  reference  vectors  x(k),  k=l,...L 
are  assumed  to  be-  independent  and  each  identically  distributed  with  zero  mean  and 
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covariance  Rsd-  The  pdf  for  the  transformed  covariance  matrix  is  also  characterized  by  the 
Wishart  distribution  [10] 


/«.(«)  = 


^  Qxp{-tr[{u^ u) 


^l/2(V-l)(V-2)p(^^  ...Y(L-N  +  2)yu 


(11) 


Now  consider  the  pdf  of  Under  conditioning,  4  is  complex  Gaussian  distributed  as 


\R, 


tz 


(12) 


n 


{N-\) 


where 

<  =  E{^^]=E{R:'VB,RTR.RTByK''^] 

Then  the  unconditional  pdf  is 


(13) 

where  (11)  and  (12)  are  inserted.  Replacing  ^^Rtz^z  with  tr[R,z  44^7  (^2)  further 

simplifies  (13).  Using  (13)  and  (9),  the  unconditional  pdf  of  the  test  statistic  is 

(14) 

cdl^s 
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Distribution  of  p 


It  turns  out  that  it  is  possible  to  find  the  distribution  of  p  in  (10)  a  different  way  which 
leads  to  an  easier  way  to  compute  quantities  like  that  computed  in  (14).  To  simplify 
things  a  bit,  let  us  ignore  rank  reduction  in  the  rest  of  this  document  so  now  we  assume 
M=N-1.  Now  consider  the  random  variable  (z  is  a  length  N-1  vector) 


P  = 


2 


z"R;'z  _  5, 
2z\  ~S, 


(15) 


Now,  from  (15),  P  is  the  ratio  of  two  random  variables.  From  [13,  p  927]  the  first,  5/.  is  a 
central  chi-squared  random  variable  with  2(L-N+2)  degrees  of  freedom  and  S2  is 
independent  of  Sj.  The  pdf  of  5y  is 

L-N+l  /■  N 

4  (-^i)  -  _  JV  +  2)  T  j 

provided  si  is  positive,  otherwise  the  pdf  is  zero.  In  the  case  without  mismatch,  from 
[7,8],  the  second  random  variable,  ^2.  is  a  central  chi-squared  random  variable  with  2(N- 
1)  degrees  of  freedom.  In  cases  with  mismatch  the  distribution  for  S2  is  more 
complicated  as  we  now  demonstrate  by  extending  some  results  from  [17].  Since  z  is 
complex  Gaussian  with  zero  mean  and  covariance  matrix 


we  can  write 


-1/2 


Pr[z“z<y]=  J 

By  factoring  Rmm=EEP  and  defining  a- 


— l^exp(-  z“R-J„z)dz 

=E'^z  we  find 
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Pr[z''  z  <  y]  =  j  ^exp(-  a''a)rfa 

a’'E'^Ea<y  ^ 


Now  by  forming  an  N-1  by  N-1  unitary  y/  =  I)  matrix  of  eigenvectors  y/  such  that 

(Ef^  E)  y/=  <P  where  0  is  diagonal  with  elements  ^ . along  the  diagonal  and 

defining  w  =  y/^  a  =  we  obtain 


f  N-\ 


Pr[z"z<>']=  J  "^exp(-w^vv)rfw  =  Pr 

w^^w<y  ^  V 


> 


Note  that  the  eigenvalues  ^>7,  02 . 0v-i  of  E^  E  are  exactly  the  same  as  those  of  Rmm 

[17].  Thus  the  distribution  of  S2  is  the  same  as  the  distribution  of 


2X0  Jw; 


where  the  WjS  are  independent  and  identically  distributed  zero-mean  and  unit-variance 
complex  Gaussian  random  variables  and  07,  02,...,  0v.7  are  the  eigenvalues  of  Rmm-  Thus  S2 
has  characteristic  function  (t  is  the  frequency  variable) 


N-l 


expUt2^^(l)j\Wj 


(17) 

where  i= V-i.  By  expanding  (17)  into  a  partial  fraction  expansion  one  finds  the  pdf  of  S2. 
First  assume  that  all  the  eigenvalues  are  distinct  so  (17)  becomes,  after  partial  fraction 
expansion 
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(18) 


N-\ 

n 


1 


1  —  2it(l). 


;=i  1  —  2it<j). 


with 


k=\,k^j 


So  the  pdf  for  S2  is  (a  mixture  of  chi-squared  random  variables  with  2  degrees  of 
freedom,  each  scaled  by  ^;) 


r 


exp 


20, 


(19) 


as  long  as  $2  is  positive  (the  pdf  is  zero  otherwise).  If  any  of  the  eigenvalues  are  repeated, 
this  is  easy  to  handle.  If  a  particular  (t>j=<p  is  repeated  r  times  this  leads  to  replacing  r 
terms  on  the  right  hand  side  of  (18)  with  (the  new  partial  fraction  expansion) 


r 

1 

1=1 


B 

_ J 

(l-2it(t>y 


and  thus  these  terms  contribute  a  weighted  sum  (weighted  by  Bj)  of  chi-squared  random 
variables  with  different  degrees  of  freedom,  each  scaled  by  or  in  equation  form  they 
contribute 


(20) 


to  (19)  as  long  as  S2  is  positive  (the  pdf  is  zero  otherwise).  Using  (19)  and  (20)  as 
appropriate  one  can  find  the  pdf  of  S2. 
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Now  returning  to  (15)  we  find 


Mp)= J"  .  ^^fs,  {p^i  )4  (■'2 


for  p  >0.  (21)  is  derived  in  the  Appendix.  For  example,  if  no  eigevalues  are  repeated  (21) 
becomes  from  (16)  and  (19) 

^  (p.,  exp(-  psj  2)^^  fc,  f  s,  V 

Jf^P)  4=0  2  2"-^"T(L-N  +  2)  I?  20.  \ 

-  2--2  r(^  _  ^  +  2)  ^  20.(/7+i/0j""' 

+1/^2))/’ + 1  /  ^2)^2 


■■2''**T(Z.-A?  +  2)^2^,(p  +  l/0j 


_  (Z,  -  W  +  2)*.p‘' 


where  we  used  the  normalization  of  a  chi-squared  pdf.  If  an  eigenvalue  (j)  is  repeated  r 
times  then  (16),  (20)  and  (21)  lead  to  a  result  with  terms  like 


r  -2^ 

Js2=o  ^  9 


(p5,)  exp(  ps2/2)-y  Bj  r  h 

2‘-«>^r(L-N+2)  p,rr{j)(i,[ii>)  H  2^ 


.a2‘“''*2«p(£_^  +  _,  +  2)  Bjp'--''*' 

L  2'-''*'r(L-W  +  2)  22T0>J(p  +  1/^)'''"'*'^ 

r  - rexpf-^(/7-Hl/0)l(;2 -1-1/0^52 

J^=02"-^"^""r(L-iV-h7  +  2)  ^ 

L  2'-"'"^"r(L-A-h2)  2'r(;)0'(/?-l-l/(/>)''"''^'^ 

_r  Bj{L-N  +  j  +  l)lp^-^^^ 

^(y_l)!(L-iV  +  l)!0^(p-Hl/0)^-"^^"^ 
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for  the  terms  corresponding  to  a  repeated  eigenvalue.  Using  (22)  and  (23)  allows  one  to 
find  a  closed  form  expression  for  the  pdf  of  P  in  any  given  case.  One  must  simply 
substitute  the  terms  in  (23)  into  (22)  in  replacement  of  the  terms  in  (22)  that  correspond 
to  repeated  eigenvalues. 

Now  we  can  use  the  distribution  of  P  to  find  the  distribution  of  p.  First  consider  the  pdf 
of  /?  =  1/P  which  is  (using  standard  transformation  of  random  variable  techniques) 

1 

fR{r)  =  fp  -  ^  (2'*) 

\r)r 

which  applies  for  positive  r  (otherwise  the  pdf  is  zero).  Then  the  pdf  of  Y=Var(d}+R  is 

fyiy)=  A {y-Var{d})  =  { - f, [ - 

\y-Var{d))  yy-Var{d}} 

(25) 

for  y  between  Var{d}  and  infinity.  Finally  p  =  1/Y  so^ 

/p(pl^^o)=/p(pl^J  =  A//-^ 

p 

fp 

for  0<p<l/Var{d},  otherwise  the  pdf  is  zero. 


^  In  this  case  the  pdf  of  p  is  the  same  under  Ho  and  Hj,  but  in  cases  considered  later  this  is  not  true.  This  is  the  reason 
for  the  notation  used  in  (26). 


p 

1-Var{d}p 
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A  Consistency  Check 

Note  that  in  the  case  without  mismatch  (Rsd  =  Rt)  we  have  ^/  =  ^2  =-••=  <I>n-i  =i  and  so 
(19)  becomes  (20)  with  r=N-l  and  Bj=B2=—=Bn-2=0,  Bn.i=1  and  so  in  this  case  S2 
becomes  a  chi-squared  random  variable  with  2(N-1)  degrees  of  freedom.  Under  this 
condition  the  random  variable  P  is  a  ratio  of  two  chi-squared  random  variables  and  from 
[8]  we  find  P  is  an  F  random  variable  with  pdf  (similar  to  (9)) 

(/.-w+i) 

iL-N  +  l){N-2)l{l  +  pY^' 

since 

^F^{L  +  l■,L-N  +  2;Q)  =  l 

Equation  (27)  applies  for  p  between  0  and  infinity  (otherwise  the  pdf  is  zero).  Now 
applying  (24),  (25),  and  (26)  we  can  use  the  transformation  to  R  and  Y  and  finally  p  for 
Var{d}=l  to  obtain  (for  0<p<l) 


4(piH„)=4(piff,)= 


{L  -  N  +  1)'{N  -2)\ 


(28) 


which  matches  the  result  in  [3]  as  expected.  As  a  second  check,  consider  the  case 
without  mismatch  again.  We  note  that  in  this  case  (23)  applies  with  (pi  =  ^  =...=  =i 
and  r=N-l  and  Bi=B2=...=Bn.2=0,  Bn-i=1.  So  we  evaluate  (23)  under  these  conditions 
and  we  obtain  (27)  and  ultimately  (28)  which  matches  the  result  in  [3]  as  expected. 

Probability  of  detection  and  false  Alarm 


The  distributions  of  p  found  in  the  last  section  lead  to  an  easier  way,  as  compared  to  (14), 
to  derive  the  detection  and  false  alarm  probabilities.  The  approach  is  similar  to  that  taken 
in  [3]  for  cases  without  mismatch.  Conditioned  on  p,  we  have  already  shown  (see 
discussion  between  (8)  and  (9))  that  our  test  statistic  has  exactly  the  same  form  as  the  one 
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in  [3,  equation  (27)]  and  so  by  the  same  arguments  given  there  we  find  that  the  false 
alarm  probability  conditioned  on  p  is 


(l  +  rp) 


L-N+l 


where  t  is  the  threshold  (3)  is  compared  to^.  It  is  interesting  that  this  is  exactly  the 
expression  for  a  scalar  CFAR  case  with  L-N+l  reference  samples  (Kelly  has  discussed 
this  in  [14]).  The  unconditional  false  alarm  probability  is 

\/Var{d}  j 

^FA  ~  J  A  \L-N+l  f  p  (P  ^  ^0 

i  (i+Tp) 

where  (26)  should  be  inserted.  The  expression  in  (30)  can  be  evaluated  by  numerical 
integration.  Under  the  detection  probability  conditioned  on  p  is  obtained  from 
[3,equation  (33)]  as 


(l+rp) 


L-N+l 


T,  (vro. 

^  (L  -  iV  + 1  -  m)\m\  1  +  'rp 


where 


as  defined  in  [3].  Using  (26)  and  (31)  gives  the  unconditional  probability  of  detection  as 


The  same  threshold  is  compared  to  (1). 
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\IVar{d} 

Pd=  (32) 

0 

Effects  of  Mismatch:  Numerical  results 

Now  we  wish  to  investigate  the  effects  of  mismatch  and  at  the  same  time  we  check  the 
correctness  of  (30)  and  (32).  For  the  case  with  mismatch  (26)  implies  that  performance  is 
sensitive  to  how  much 

Var{d]=v’'RTR.RT^  (33) 

differs  from  unity.  Further,  (22)  and  (23)  indicate  that  performance  is  also  sensitive  to 
how  much  the  eigenvalues  of 

p..=b.RTR.RTb'!  (34) 

differ  from  unity.  Varfd}  and  the  eigenvalues  of  Rmm  are  the  only  parameters  which 
effect  performance'^  so  we  present  results  for  various  values  of  these  parameters.  We 
compare  calculated  results  using  (30)  and  (32)  to  Monte  carlo  simulations  using  50,000 
trials. 


Consider  the  specific  case  of  N-=4,  L=8  and  a  probability  of  fedse  alarm  of  0.01.  Similar 
results  are  obtained  in  other  cases.  First  consider  the  case  with  the  eigenvalues  of  Rmm 
fixed  at  C=/^/,  ].  The  results  in  Fig.  1  show  that  the  probability  of  detection 

generally  increases  as  Var{d}  decreases  with  C  fixed.  This  is  reasonable  since  1/Var{d} 
measures  the  SNR  of  d  and  any  signal  energy  will  end  up  in  d.  Note  the  excellent 
agreement  between  the  curves  produced  by  (30)  and  (32)  and  the  “points”  produced  by 
the  Monte  carlo  simulations.  In  Fig.  2,  Var{d}  is  fixed  and  C  is  varied.  We  see  that 


^  Of  course  performance  also  depends  on  N  and  L  as  indicated  in  (22)  and  {23)»  leading  to  a  dependence  in  p  similar  to 
that  in  (27). 
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increasing  any  eigenvalue,  is  similar  to  increasing  Var{d}.  This  is  reasonable  since 
each  eigenvalue  can  be  thought  of  as  the  noise-plus-clutter  power  from  unwhitened  noise- 
plus-clutter  in  a  particular  one-dimensional  subspace.  Clearly  increasing  noise-plus- 
clutter  power  should  lead  to  a  decrease  in  performance.  Again,  there  is  excellent 
agreement  between  the  calculated  and  simulated  results.  In  the  results  in  Fig.  1  and  Fig. 
2  the  threshold  of  the  test  is  changed  each  time  Var{d}  or  C  is  changed  so  that  the 
probability  of  false  alarm  is  fixed  at  0.01.  This  provides  an  easy  to  interpret  result. 

In  practice,  both  probability  of  detection  and  probability  of  false  alarm  will  be  altered  due 
to  mismatch  so  now  we  consider  changes  in  false  alarm  probability.  Here  we  fix  the 
threshold  at  the  value  needed  to  achieve  the  correct  probability  of  false  alarm,  0.01,  for 
the  case  without  mismatch,  Var{d}=l  and  C=[l,l,l].  Fig.  3  gives  results  for  fixed  C  and 
Fig.  4  gives  results  for  fixed  Var{d}.  The  results  indicate  that  increasing  Var{d}  or  any 
of  the  components  of  C  leads  to  an  increase  in  the  probability  of  false  alarm.  This  can  be 
expected  based  on  the  interpretations  of  Var{d}  and  C  in  the  last  paragraph.  Increasing 
either  of  these  is  like  increasing  the  effective  noise-plus-clutter  power. 

Now  we  consider  an  airborne  radar  example  where  ground  clutter  produces  a  clutter  ridge 
in  the  two  dimensional  spatial-Doppler  power  spectral  density  (psd).  For  simpilicity,  we 
consider  /?,  and  R^d  modeled  by  the  simple  model  described  in  [11].  This  model  places 
Gaussian  humps  along  a  line  in  the  angle-Doppler  plane  to  model  the  clutter  ridge  as 
shown  in  the  psd  in  Fig.  5.  The  equation  for  this  psd  is 


w,./,)=y 


<T 


c,d 


exp 


r 


2<t 


ft  4 


2(J 


fs4 


which  is  a  function  of  doppler  frequency  and  spatial  frequency  .  This  psd  consists 

of  L  Gaussian-shaped  humps,  the  d  th  of  which  is  centered  at  » /, )  =  (/«  d->fcsd) 

2 

and  has  amplitude  controlled  by  (7^^  and  a  spread  in  angle  and  doppler  controlled  by 
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2  2 

and  (7^  ^ .  The  choices  of  the  parameters  used  in  the  above  equation  are  given  in 

Table  1.  The  extra  hump  in  Fig.  5  is  inserted  just  in  either  the  psd  for  the  cell-under-test 
(corresponds  to  altering  R,  [11])  or  in  the  psd  for  the  reference  data  (corresponds  to 
altering  Rsd  [11]).  This  extra  hump  models  mismatch,  possibly  due  to  some  discrete 
scatterers.  The  target  is  located  near  the  extra  hump  to  have  maximum  effect.  The 
performance  results  are  shown  in  Fig.  6  for  several  different  choices  of  target  location 
and  extra  hump  location.  For  each  different  choice  for  the  location  of  the  extra  hump,  the 
threshold  of  the  test  is  changed  to  maintain  the  false  alarm  probability  of  0.01.  In  all 
cases,  if  the  hump  is  added  in  the  reference  data  performance  improves,  but  if  the  hump  is 
added  to  the  cell-under-test  data  performance  degrades.  This  simple  example  illustrates 
the  use  of  our  equations  in  an  airborne  radar  application. 

One  can  explain  the  results  in  Fig.  1  through  Fig  .  6  by  considering  the  changes  in  the  pdf 
of  p  caused  by  changes  in  var{d}  or  C.  Any  of  the  decreases  in  Var{d}  or  ())j  causes  the 
mass  in  the  pdf  of  p  to  move  towards  larger  values  of  p.  Due  to  the  decreasing  nature  of 
the  function  multiplying  the  pdf  of  p  in  the  integrand  of  (30),  this  causes  a  decrease  in  the 
probability  of  false  alarm.  It  is  possible  to  apply  similar  analysis  to  explain  changes  in 
probability  of  detection  also. 

Extension  to  a  more  general  class  of  Tests 


Consider  the  class  of  tests  considered  in  [12]  which  compare  Aamf  to 

Y[1  +  COC^R~^x\  (35) 

If  a=0  this  is  the  test  we  considered  previously  with  t=)'  If  a=i  then  this  becomes  the 
Generalized  Likelihood  Ratio  (GLR)  test  [4].  In  [12]  a  is  called  a  sensitivity  factor.  The 
only  difference  between  the  above  formula  and  (2)  in  [12]  is  due  to  the  missing  1/L  in  our 
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(2)  (compare  to  (3)  in  [12]).  Applying  the  same  transformations  as  before  yields  that  this 
is  equivalent  to  comparing  the  simplified  test  statistic  in  (3)  to  the  threshold 


7 

l-ya 


(l-a)  +  ^ 

P. 


(36) 


This  means  we  can  find  the  probability  of  false  alarm  and  the  probability  of  detection  for 
this  more  general  class  of  tests  by  using  (29)  through  (32)  with  t  replaced  by  the  quantity 
in  (36). 


It  is  of  interest  to  extend  the  results  to  other  STAP  algorithms.  One  extension  is 
considered  in  the  next  section, 


Steering  and  Covariance  Mismatch 


Now  consider  the  case  where  both  the  steering  vector  and  the  covariance  matrix  have 
mismatch.  This  case  provides  performance  for  ADPCA  for  cases  with  covariance  matrix 
mismatch.  Consider  the  test  in  (1)  when  the  steering  vector  is  q,  but  were  the  true  signal 
vector  that  is  observed  is  the  vector  Ks.  As  in  the  previous  case,  first  apply  a  coordinate 
transform,  which  consists  of  multiplication  by  to  the  observed  vector  from  the 

cell-under  test,  the  reference  data  vectors  and  to  the  signal  and  steering  vectors.  Next, 
normalize  the  transformed  signal  vector  so  that  it  is  again  a  unit  vector.  Call  the 
transformed  unit  signal  vector  v.  So  that  signal-to-noise  ratio  is  maintained,  the 
transformed  signal  is  taken  to  be  [(s^  Rsd  kJ  v.  Similarly,  call  the  normalized 
transformed  steering  vector  qu.  Next,  another  transformation  of  coordinates,  which 
consists  of  multiplication  by  (  qu,Bq^)^,  is  made  to  the  cell-under-test  vector,  the 
reference  data  vectors  and  to  the  signal  and  steering  vectors.  Here  Bq  is  a  matrix  whose 
N-1  rows  consist  of  a  set  of  N  dimensional  vectors  that  span  the  space  orthogonal  to  qu. 
The  matrix  (  qu,B^}^  is  taken  to  be  unitary.  Denote  the  transformed  vector  for  the  cell- 
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under-test  by  y  =  ( d,  f-  Here  d  is  a.  scalar  that  describes  the  component  of  the  vector 
Rsd  '^x  which  lies  in  the  direction  of  qu-  On  the  other  hand,  z  is  an  N-1  dimensional  vector 
describing  the  component  of  Rsd  ^^x  which  is  orthogonal  to  qu-  For  the  transformed 
reference  vectors  employ  the  notation  y(k)  =  (  d(k),  z(kf  f.  Under  the  assumptions 
outlined,  we  again  obtain  (3).  Further,  the  previous  observations  about  (3)  still  hold, 
except  that,  under /f/,  the  expected  value  of  yr  given  z,  z(l),  z(2),...,z(L)  is  now 


(37) 


and  so 


Thus  the  result  in  (38)  is  similar  to  the  previous  result  (see  the  distribution  of  Aamf  just 
prior  to  (9)),  but  with  a  different  effective  value  of  fi.  The  results  like  (9)  and  (15)  still 
hold  (with  P  set  as  suggested  by  (38))  and  in  fact  (26)  still  gives  the  correct  distribution 
of  p  under  Ho  if  we  define 


However  (26)  does  not  give  the  correct  distribution  of  p  under  Hi.  Now,  since  q  and  s 
don’t  point  in  the  same  direction,  only  part  of  the  signal  goes  into  d,  while  the  other  part 
goes  into  z.  Thus,  under  Hi  we  find  S2  =  2  z^z  contains  signal  which  alters  the  results 
given  by  (19)  and  (20)  and  thus  (22)  and  (23).  This  is  the  key  difference  in  this  case. 
The  signal  in  z  also  changes  the  effective  SSNR.  As  we  see  in  (38)  the  effective  value  of 
SSNR  is  reduced  from  IjSl^to 


1C  k  =\Pfcos\e) 

which  might  be  called  the  SSNR  in  the  subspace  spanned  by  the  steering  vector.  Here 
cos(Q)  is  defined  by  this  equation  using  the  previous  definition  of  \p\^=\KS^(s^ s).  For 
simplicity  we  take 


where  an  unimportant  possible  complex  multiple  has  been  ommitted  [3].  Now  the  part  of 
lj3P  which  was  lost  went  into  z  and  is 

kf  =  \E{z  \hJ=  \pf  sin"(0)  = 

which  could  be  called  the  SSNR  in  the  subspace  orthogonal  to  the  steering  vector. 

The  signal  entering  z  changes  the  distribution  of  S2.  Thus,  we  now  consider  the 
distribution  of  S2  (the  distribution  of  Sy  is  still  given  by  (16)).  Again  we  extend  some 
results  in  [17].  Under  Hj  z  is  Gaussian  with  mean  =BqRsd^Ks  and  covariance  matrix 
Rmm-  Thus,  we  can  write 

Pr[z''z<>']=  J  ;3r4^exp(-(z-/iJ''/f;Uz-A)Vz 

By  factoring  Rmm=^^  and  defining  a=E'^z  and  iia=  we  find 

Pr[z''z<3']=  J  -4rexp(-(a-/ij"(a-ytj)tffl 
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Now  by  forming  an  N-1  by  N-1  unitary  (y/^y/  =  I)  matrix  of  eigenvectors  y/^such  that 

y/^  E)yf=  0  where  0is  diagonal  with  elements  0/,  ^ . <^^-1  along  the  diagonal  and 

defining  u=y^  a  and  od=  y/^ iia  =  (C0j„..,CDN.]f  vfe  obtain 


Pr[z"z<y]=  j  _Lexp(-(M-ffl)''(M-<o))tf, 

u»t>u<y  ^ 


u 


J  -IrexpC-  w"  w)rfw  =  Prf  +  (o. 

rh  ( W  ■»»  71  y  y=l 


< 


J 


where  w  =  u-co.  Again,  we  note  that  the  eigenvalues  <p2 . ^/v./  of  E  are  exactly 

the  same  as  those  of  Rum  [17]. 

Thus  the  distribution  S2  =  2  z^z  is  the  same  as  the  distribution  of 


(39) 


where  as  before  the  WjS  are  independent  and  identically  distributed  zero-mean  and  unit 
variance  complex  Gaussian  random  variables.  In  (39)  (CDi,...,  (%-//  E'^BvRsd^^Ks. 
Also,  as  before  0;,  are  the  eigenvalues  of  Rmm-  Now  the  characteristic  function 

for  S2  must  be,  using  (39) 
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ji[  j  j  exp{lit(j)j  [(x  +  6)^.  y  +  (y  +  (0,j  J  ])—  exp(-  {x^  +  ))dxdy 

;^i  7C 

J-i  x=-oo  y=-c« 


£_^exp(-  (>>  +  fiij  y  (l  -  2it<pj  )}ly 

exp(^it(p\co^f)sxpl^l  +  )) 


r 


=  J|(l-2iY0^.]r‘exp 


y=i 


V 


2it(l)j  (O. 
1  —  2it^. 


|2  \ 


(40) 


where  (Oj  =  (Orj  +  icoij  and  fXj  =  fiRj  +  ifiij  are  suitably  chosen  to  make  (40)  true  (pick  the 
cross  terms  in  the  quadratic  in  the  exponential  to  match).  Now  note  from  (40),  that  if  all 
=  1  and  (Oj  =  0  (this  occurs  without  mismatch,  where  Rsd=Rt  and  s=q)  then  (40)  implies 
that  52  is  central  chi-squared  with  2(N-1)  degrees  of  freedom  as  suggested  by  the  known 
result  for  the  case  without  mismatch.  If  there  is  mismatch  in  covariance,  but  no 
mismatch  in  Steering  vector,  then  (c^  =  0  and  (40)  matches  (17).  Suppose  there  is  no 
mismatch  of  covariance  matrix  R^d—Rt  but  that  there  is  a  mismatch  in  steering  vector  so 
syiq.  Then  from  (40),  S2  has  a  noncentral  chi-squared  pdf  with  2(N-1)  degrees  of 
freedom,  which  is  denoted  by  ^2^iv.//l<ol^).  Thus  R=l/P  has  an  F  distribution  [8] 
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.  s  L!r'""^exp(-b^  )  f  i 

(r  1 //,  )=-^ - — - r-: - ^1^1  ^  +  1;A^— 1;^:; - 

jRiH,\  iJ  (^i _js/  +  i)(jv- 2)1(1 +  ry  1-rJ 

The  pdf  of  p  can  be  found  using  (25),  (26)  and  the  series  expansion  for  /Fy  in  [8]  to  be 


f  (jin)  g^p(-KI) f ^  I  i-jv,  1- 


=  exp(-KI)t  IbT  ^ - (i±]±£!> - p) 

^  '  '’timf}  ’  (Z.  +  l-iV>(Ar  +  m-2) 

which  checks  with  (30)  from  [3]  which  is  appropriate  for  this  case. 


To  invert  (40)  to  obtain  a  pdf  (in  the  general  case)  it  is  useful  [17]  to  expand  (40)  in  terms 
of  6=  1/(1  -  2itQ  by  rewriting  (40)  as 


E<  exp  Wj  +  cOj 


exp 


1-e 


2k»: 


2  J  M  c  P " 

1-0  1--^ 

d). 

V 


Hi- 1-^ 


By  canceling  0^'^  from  both  sides  of  (41),  the  coefficients  of  the  expansion  can  be  found 
by  using  the  standard  formula  for  Taylor  series.  They  can  be  expressed  as 
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«o=exp-XK 

e  = - ^  B  for  r>0 


2r^ 


A^-i  ^  ^  Y  N-ij  26)  .  ^  f  ^ 

B,  =  l2  1-f  +rX-f^  1-f 

M  <Pj  >.  (Pj  ,  <I>J 


Now  note  that  (41)  says  we  can  represent  the  distribution  of  S2  as 

/  \N-2+j 


fs  7  '■ - r~GXp 

Js,\2J  ^  >  2^-'*'r(N -1+ j)g 


which  is  a  mixture  of  central  chi-squared  pdfs.  In  (45),  we  set  ^=1  if  all  the  (j)j  >1. 
Otherwise  we  set  ^  =  minj  <t)j.  Now  use  (16),  (21)  and  (45)  to  find  that  the  pdf  of  p  is 


in'l-f  ? 

ipS2. 

\L-N+l  ( 

)  exp( 

2^- 

■""TfL- 

.N-2+j 


OQ 


fs  '2''-'*'r{N-l  +  j)g 

which  simplifies  to 


rr-exp  - 
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/f(p)=Z«; 

7=0 


p‘-h»2^*‘*'T{L+  ;  +  1) 


2“"'  T{L-N  +  2)T{N  - 1  +  jY  p  +  ^ 


\L+;+l 


X 


f 

J57=0 


oo 


p+ 


exp 


-  i7  + 


V.v- 


aAj 


p‘-^*'T{L+ j +  1) 


!'  1 
P  +  ^ 


2=0 


r(L-Af  +  2)r(Af-l  +  7)^p  +  ^ 


nL+7+1 


.^-1+7 


(47) 


for  p  >0.  Using  (47)  in  (26)  provides  the  desired  distribution  for  p  under  Hj.  Then  (29), 
(30),  (31)  and  (32)  provide  the  probability  of  false  alarm  (using  the  previously  found  pdf 
for  p  under  Ho)  and  detection. 


Now  note  that  if  =  02  =•••=  <Pn-i  =1  (which  occurs  if  Rsd  =  Rt,  thus  no  mismatch) 

then  eo  =  1  and  ej  =  0,j  ^0  from  (42),  (43),  and  (44).  In  this  case  (45)  becomes  a  2(N-1) 
central  chi-squared  pdf  as  required  and  (47)  matches  with  (27). 

Steering  and  Covariance  Mismatch: 

An  Explanation  for  the  good  performance  ofADPCA 

We  have  performed  numerical  studies  for  cases  with  steering  vector  mismatch  similar  to 
those  in  Fig.  1  through  6.  These  results  illustrate  that  performance  varies  with  Var{d} 
and  <j)j,  j=l,2,...,N-l  in  a  manner  similar  to  that  shown  in  Fig.  1  through  6.  Also 
plays  a  role  similar  to  that  played  by  1)81^  (SSNR)  in  Fig.  1  through  6.  Finally,  tu  is  a  new 
parameter  for  the  case  with  steering  vector  mismatch,  but  it  is  easy  to  understand  the  role 
of  a>.  An  increase  in  any  component  of  (o  generally  leads  to  decreases  in  probability  of 
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detection  and  no  change  to  probability  of  false  alarm.  Since  co  is  a  transformation  of  the 
signal  that  is  orthogonal  to  the  steering  vector  this  is  reasonable. 

Probably,  the  most  interesting  aspect  of  the  results  is  that  steering  vector  mismatch  can 
sometimes  offset  covariance  matrix  mismatch.  In  fact,  this  appears  to  explain  the  good 
performance  of  ADPCA  which  was  observed  in  pervious  studies.  Recall  that  one  key 
aspect  of  ADPCA  is  that  it  does  not  use  the  optimum  steering  vector  for  the  case  where 
the  covariance  matrix  is  exactly  known.  Instead,  it  uses  a  steering  vector  that  can  whiten 
clutter  that  was  not  completely  whitened  after  multiplication  by  the  inverse  of  the 
estimated  covariance  matrix.  Some  examples  of  cases  where  steering  vector  mismatch 
offsets  covariance  matrix  mismatch  are  provided  in  Fig.  7.  In  each  case  we  see  the 
probability  of  detection  curve  for  the  case  of  both  steering  and  covariance  matrix 
mismatch  is  above  the  curve  for  the  case  with  just  covariance  matrix  mismatch.  We  note 
that  this  does  not  always  occur.  In  Fig.  8,  we  show  cases  where  the  probability  of 
detection  curve  for  the  case  of  both  steering  and  covariance  matrix  mismatch  is  below  the 
curve  for  the  case  with  just  covariance  matrix  mismatch. 

In  fact,  it  is  easy  to  explain  one  case  where  steering  vector  mismatch  can  sometimes 
offset  a  covariance  matrix  mismatch.  Consider  the  case  where  L  is  very  large  so  that  Re 
is  very  close  to  i?,  in  the  case  without  covariance  matrix  mismatch  and  Re  is  very  close  to 
Rsd  in  the  case  with  mismatch.  Now  recall  that  Re  is  used  in  the  test  statistic  in  (1)  and 
that  if  L  is  very  large  the  denominator  of  (1)  is  close  to  a  constant  which  can  be 
considered  part  of  the  threshold  that  is  chosen  to  set  the  required  false  alarm  probability. 
Now  it  is  easy  to  imagine  that  one  can  choose  a  mismatched  q  that  can  make  Rsd^ 
closer  to  the  optimum  value  of  Rf^  than  Rsd'^  is.  In  fact,  this  is  exactly  the  case 
shown  in  Fig.  6. 
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Conclusion 


We  would  like  to  continue  studying  mismatch  in  order  to  develop  STAP  schemes  that  are 
robust  to  mismatches.  We  believe  this  is  an  important  practical  issue  that  deserves 
attention  and  that  significant  advanced  can  be  made  through  such  research.  We  believe 
this  is  clearly  demonstrated  by  the  results  given  in  this  report. 
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Appendix 


SimpliHcation  of  Aamf 

Apply  the  transformations  and  normalizations  as  required  to  (1).  Then  use  the  Frobenius 
relation  [4]  in  (1)  to  obtain 


hi  -  ri rj'd-  r “  Rj'  (a]  -  Rj' )"'  J 


-r5R; 


Simplification  of  this  (multiply  numerator  and  denominator  by  the  scalar  equal  to  the 
square  of  the  inverse  of  the  denominator)  gives  (3). 

Rank  Of  C 


Defining 

B  =  (fc(l),fe(2), I.))" 

and  assuming  L>M  and  that  B  is  full  column  rank  (this  occurs  with  probability  1)  gives 

=  (u^b) 

Now  noting  that  the  column  rank  of  u  is  M,  implies  that  the  rank  of  all  of  the  three 
matrices  in  the  product  above  must  be  M  and  so  the  product  must  have  rank  M  [8, 

Appendix].  Thus  multiplication  by  Z^R~^Z  projects  onto  an  M  dimensional  space. 

Similarly  /  —  Z^^”'Z  (7  is  an  L  by  L  identity  matrix)  projects  onto  an  L-M  dimensional 
space. 


Idempotence 


By  direct  calculation 
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(z"^;'z)(z"^;'z)  =  (z"^;'z) 

Variance  of  yr 

Using  (8),  the  independence  of  d,d(l d(L),  and  conditioning  on  z,  z(l z(L)  gives 

(  I 

Var{y;}=Var{d}+E\  '^d(k)z(k)“ R:'z 

=  Var{d}+  =  Var{d]+  E^(kf}z“  R;' z 

since  d(l),  .....  d(L)  are  iid.  Also 

and 

E{<i(t)"}  =  l 
Derivation  of  (21) 

f  \  oo  / (p+dp)s2  ^ 

Pt\p<-^<  p  +  dp  =f^{p)dp  =  j  ds^dp 

^2  )  0\ps2  ) 

OO 

=  J {vh^h)hds^dp  =  J 4^  {ps,)f^  {s;)s,ds,dp 

0  n 
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Probabriily  of  Detection  Probability  of  Detection  ProbabiNty  of  Detection  ProbaWlity^  D^ection 


SSNRIdB) 


Fig.  1.  Probability  of  detection  variation  with  Var{d}. 
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Probability  of  Detection  Probabilily  of  Detection 


SSNR(dB)  SSNR(dB) 


Fig.  2,  Probability  of  detection  variation  with  C=[  ,  (j)^,  (j)^]. 
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15 


Fig.  4.  Probability  of  false  alarm  variation  with  C. 
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Fig.  5.  Power  spectral  density  of  a  simple  model  [17] 
for  ground  clutter  and  mismatch  (extra  hump). 


Table  1.  Parameters  of  the  power  spectral  density  model  in  Fig. 


d 

f rt,rf 

f ts4 

1 

-0.35 

-0.35 

0.03 

0.03 

8.00 

2 

-0.20 

-0.20 

0.03 

0.03 

8.00 

3 

0.00 

0.00 

0.03 

0.03 

9.98 

4 

0.20 

0.20 

0.03 

0.03 

8.00 

5 

0.35 

0.35 

0.03 

0.03 

8.00 

extra  hump 

0.38 

0.18 

0.03 

0.03 

7.00 
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Fig.  6.  Performance  variation  for  an  airborne  radar  example 
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ProUit/ of  detection  Probilit/ of  detection  Probility  of  detection 


Fig.  7.  Probability  of  detection  variation  with  steering  mismatch  (offset). 


Fig.  8.  Probability  of  detection  variation  with  steering  mismatch  (not  offset). 
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